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Abstract. Some of the most 'active' galaxies in the Universe are obscured by large 
quantities of dust and emit a substantial fraction of their bolometric luminosity 
in the infrared. Observations of these infrared luminous galaxies with the Infrared 
Space Observatory {ISO) have provided a relatively unabsorbed view to the sources 
fuelling this active emission. The improved sensitivity, spatial resolution and spectro- 
scopic capability of ISO over its predecessor Infrared Astronomical Satellite {IRAS), 
has enabled significant advances in the understanding of the infrared properties of 
active galaxies. ISO surveyed a wide range of active galaxies which, in the context 
of this review, includes those powered by intense bursts of star-formation as well as 
those containing a dominant active galactic nucleus (AGN). Mid infrared imaging 
resolved for the first time the dust enshrouded nuclei in many nearby galaxies, while 
a new era in infrared spectroscopy was opened by probing a wealth of atomic, ionic 
and molecular lines as well as broad band features in the mid and far infrared. This 
was particularly useful since it resulted in the understanding of the power produc- 
tion, excitation and fuelling mechanisms in the nuclei of active galaxies including 
the intriguing but so far elusive ultraluminous infrared galaxies. Detailed studies 
of various classes of AGN and quasars greatly improved our understanding of the 
unification scenario. Far-infrared imaging and photometry also revealed the presence 
of a new very cold dust component in galaxies and furthered our knowledge of the 
far-infrared properties of faint starbursts, ULIGs and quasars. We summarise almost 
nine years of key results based upon ISO data spanning the full range of luminosity 
and type of active galaxies. 
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1. Introduction 

"Activity" is a common feature of galaxies tiiat emit a large fraction of 
their energy in the infrared. Existing primarily in a state of quiescence, 
during certain periods of their evolution galaxies experience short-lived 
phases of extreme star formation and/or black hole activity. Prom early 
ground-based measurements (e.g. Low and Kleinmann, 1968; Klein- 
mann and Low 1970b, 1970a) to surveys performed by the Infrared 
Astronomical Satellite {IRAS, Beichman et al., 1988), many previously 
known active galaxies [e.g. optical- or radio-selected quasars, active 
galactic nuclei (AGN), and starbursts] were found to be as or more 
powerful in the infrared than in the visible. In addition, a new popu- 
lation of optically faint, infrared luminous galaxies (luminous infrared 
galaxies (LIGs, with L//? > lO^^L^) were discovered that emit most of 
their bolometric luminosity in the infrared (Houck et al., 1984). AGN 
and massive bursts of star formation are the only known mechanisms 
that are capable of generating such luminosities. However, without 
spectroscopic information of the obscured components, classifying the 
underlying power source has proved to be difficult. Furthermore, quan- 
titatively assessing the contribution to the total power from the two 
possible mechanisms is non-trivial. It is known that dust is formed 
during the late stages of stellar evolution and it is destroyed by intense 
radiation fields and shocks. Consequently, a detailed understanding of 
the physical mechanisms responsible for the observed infrared emission 
is an essential first step into revealing the nature of the underlying 
exciting source. As a result, since the most 'active' regions in our 
Universe appear to be enshrouded in large quantities of gas and dust, 
addressing the aforementioned points would have serious consequences 
in estimating star formation and black hole activity in the Universe. 

In a cosmological context, both the reported increase in the star 
formation density of the Universe for 0.1 ^ z ^ 2 (e.g. Madau et al., 
1996; Barger et al., 2000) and the high frequency of starbursts hosted 
in galaxies with disturbed morphologies or in interacting/merging sys- 
tems, imply that starbursts play an important role in galaxy forma- 
tion and evolution. In addition, black hole activity could be buried in 
any luminous infrared system. Ultraluminous infrared galaxies (ULIGs, 
LiR> 10^2^0) have been proposed to be an evolutionary stage in the 

* Based on observations with ISO, an ESA project with instruments funded by 
ESA Member States (especially the PI countries: France, Germany, the Netherlands 
and the United Kingdom) and with the participation of IS AS and NASA. 
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life of a quasar and the local analogues of the z>2 (sub-)mm popula- 
tion. Quantifying the amount of obscured activity at earlier times and 
its contribution to the infrared and X-ray backgrounds remain open 
issues requiring a sound understanding of the properties of local active 
galaxies. 

Following on from the coarse, photometric legacy of IRAS, ISO 
provided the first means to investigate the physical conditions of local 
active galaxies in significant detail at wavelengths where the obscured, 
and often most active, regions could be probed. Equipped with imaging, 
photometric and spectroscopic capabilities, ISO surveyed a wide range 
of emission properties from broad-band spectral energy distributions 
(SEDs), through imaging, to detailed spectral analysis. We discuss the 
infrared emission properties of obscured active galaxies that ISO sur- 
veyed during its lifetime. In the context of this chapter, 'active' refers to 
both star formation and black hole activity and encompasses the range 
of non- normal galaxies: interacting/merging galaxies; starbursts; radio 
galaxies; AGN; quasars; low ionisation nebular emission regions (LIN- 
ERs); and the full suite of LIGs - ULIGs and hyperluminous infrared 
galaxies (HyLIGs > lO^^^^) 



2. Activity manifest in the infrared 



High energy UV (and visible) photons emanating from active regions 
can heat or excite environmental dust, located around active sites 
or distributed throughout the interstellar medium (ISM), resulting in 
re-radiation in the infrared. Grains ranging from the very small (ra- 
dius a^lOnm) to the large (a~30/im) contribute to emission over the 
2.5-200;um range. The type, size and distribution of dust grains to- 
gether with the incident interstellar radiation field (ISRF) shape the 
observed spectral features and overall form of the SED which arise 
from both thermal and non-thermal processes. At long infrared wave- 
lengths (A ^25/im), emission predominantly arises from grains that 
re-radiate in thermodynamic equilibrium. In young radio galaxies, ther- 
mal bremsstrahlung may also provide a contribution to the infrared 
continuum. Emission related to transient rather than steady-state heat- 
ing of dust grain complexes can dominate the SEDs of some active 
galaxies giving rise to continua and features over 3-18/im. Non-thermal 
components include the featureless far-infrared tail (increases with in- 
creasing wavelength) of the synchrotron radiation spectrum in strong 
radio galaxies. 
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2.1. Infrared continua: thermal and non-thermal 

The mid-infrared spectra of galaxies display two continua: a flat con- 
tinuum for A ^5/im and, in some galaxies, an unrelated continuum for 
A ^11/Lim (Helou et al., 2000). These continua are not produced by 
grains heated to thermal equilibrium by the radiation field but arise 
from the stochastic heating of small dust grains (amorphous silicates 
and graphites, a~0.01-0.1//m) which are transiently heated by the ab- 
sorption of a single photon (Sellgren, 1984; Beichman, 1987; Boulanger 
et al., 1998b; Draine and Li, 2001; Drainc, 2003). A single photon with 
energy of only a few electron volts (A ^0.4;um), is sufficient for the onset 
of infrared emission (Uchida et al., 1998). 

The steeply rising (near-thermal) continuum longward of ll/xm, 
which is strong in some active galaxies (but weak in quiescent galax- 
ies), is thought to originate from the excess transient heating of very 
small fluctuating dust grains (VSG, a< lOnm). This component ap- 
pears to be a characteristic of intensely star forming regions (Desert 
et al., 1990; Verstraete et al., 1996; Cesarsky et al., 1996b; Laureijs 
et al., 1996). Under this transient heating regime, the mid-infrared 
continuum is directly proportional to the underlying radiation field 
over several orders of magnitude and varies with the conditions of the 
HII regions from which it originates (Boulanger et al., 1998a; Lau- 
rent et al., 2000; Forster Schreiber et al., 2004). For example, in the 
close vicinity of the dense starburst knots in the Antennae, this con- 
tinuum is shifted to shorter mid-infrared wavelengths implying higher 
temperatures (Vigroux et al., 1996). 

The 3/Ltm< A <5/Ltm continuum is likely to arise from a feature- 
less fluctuating component (Helou et al., 2000). Although often weak 
in active galaxies, the direct stellar continuum can also contribute 
(A <5//m) and should be considered in careful spectral decomposition 
analysis. The strength of this component increases with decreasing 
wavelength and thus the transition between direct stellar to interstellar 
dust emission occurs within the near- to low- mid- infrared (Boselli et al., 
1997; BoseUi et al., 1998; Alonso-Herrero et al., 2001). The mid-infrared 
also includes the transition from the transient to the steady-state heat- 
ing regime. The wavelength at which this transition occurs is deter- 
mined by the incident radiation field above which grains of a certain 
size no longer suffer large temperature changes. Dust grains are heated 
by the general ISRF and attain a 'steady-state' where emission and 
absorption reach an equilibrium. The grain temperature at which this 
equilibrium occurs characterises the emission giving rise to a thermal 
black-body continuum. Qualitatively, the mid- to far-infrared continua 
of active galaxies are often approximated by a series of superposed 
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black-bodies, each originating from a body of dust heated to a char- 
acteristic temperature (Haas et al., 1998a; Ivison et al., 1998; Klaas 
et al., 2001; Bendo et al., 2003; Stevens et al., 2005). This approx- 
imation provides only an indication of typical temperatures and is 
not an accurate description of the complex heating of a variety of 
dust grains. More accurate descriptions are achieved by considering a 
multi-grain dust model and complex radiation fields using, for example, 
semi-empirical models (e.g. Dale et al., 2001a) or radiative transfer 
theory (e.g. Acosta-Pulido et al., 1996; Silva et al., 1998; Alexander 
et al., 1999a; Sicbcnmorgen et al., 1999b; Efstathiou et al., 2000; Verma 
et al., 2002; Farrah et al., 2003; Freudling et al., 2003; Gonzalez- Alfonso 
et al., 2004; Siebenmorgen et al., 2004a). For sources that are spatially 
unresolved by the ISO instruments, SED decomposition (using one of 
the methods above) is the only means to investigate the underlying 
physical processes. 

Some AGN-hosting galaxies display a relatively broad and flat (in 
uF^) hot continuum that dominates their mid-infrared SEDs. This 
emerges from grains in the putative torus that are heated up to the 
grain sublimation temperature {Tgub ~1500K) by the strong radiation 
field of the accretion disk of the central AGN. The 'warm' IRAS galaxy 
criterion S60/S25 > 0.2) selects galaxies displaying this component 
which are mostly identified with AGN (de Grijp et al., 1985). 

The extended wavelength range of ISO over IRAS enabled (a) the 
precise determination of the turnover in the far-infrared SEDs of ac- 
tive galaxies and (b) the identification of two dust components that 
comprise the far-infrared emission that are both associated with large 
grains in thermal equilibrium. The first (or 'cold') component repre- 
sents dust surrounding star forming regions and thus its strength is 
directly related to the star formation rate (however see Sect. 4.5). 
While the SEDs of normal quiescent galaxies peak at ~150/xm, the 
enhanced heating in actively star forming regions produces a strong 
far-infrared peak at ~60-100/xm that dominates the infrared SED and 
corresponds to dust heated to 30-60K (e.g. Calzetti et al., 2000; Dale 
et al., 2001a; Spinoglio et al., 2002). This component also enhances 
the mid-infrared continuum (e.g. see Figure 1 in Sanders and Mirabel, 
1996). In AGN, the fractional contribution of this 'cold' far-infrared 
component can be weaker in comparison to the warmer mid-infrared 
component (A ~ 30/um) than in starbursts, reflective of the intense 
heating power of the AGN. The second ('very cold' or 'cirrus') compo- 
nent is dust that is spatially more extended and is associated with the 
HI disk or the molecular halo (Trewhella et al., 2000; Radovich et al., 
2001; Stickel et al., 2004). Located far and, in some galaxies, shielded 
(Calzetti et al., 2000) from the intense emission from the central engine. 
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the extended component is heated by the ambient interstellar radiation 
field to temperatures of 10-20K, (Odcnwald ct al., 1998; Alton et al., 
1998; Davics ct al., 1999; Haas ct al., 1998b; Haas, 1998; TrewhcUa 
et al., 2000). These temperatures are consistent with the theoretically 
predicted heating/ cooling of dust grains in thermal equilibrium by a dif- 
fuse radiation field. Both profile measurements for resolved sources and 
far-infrared-sub-mm SED analysis for more distant galaxies reveal that 
a very cold component is required to explain all the emission beyond 
lOOjLtm (Rodriguez Espinosa et al., 1996; Radovich et al., 1999; Sieben- 
morgen et al., 1999a; Dominguc ct al., 1999; Calzetti et al., 2000; Haas 
et al., 2000a; Klaas et al., 2001; Perez Garcia and Rodriguez Espinosa, 
2001; Haas et al., 2003). This component, that IRAS was insensitive to, 
may be a significant contributor to the far-infrared luminosity in some 
active galaxies (e.g. ^ 60% in local starbursts Calzetti et al., 2000). 

2.2. Features 

The advent of ISO spectroscopy enabled the first sensitive and high 
resolution mid- to far-infrared spectroscopic measurements of active 
galaxies providing a better understanding of the physical conditions 
within visually obscured regions of systems such as Circinus (Moor- 
wood et al., 1996; Sturm et al., 2000), the starburst M82 (Sturm et al., 
2000) and the Seyfert 2 NGC1068 (Lutz et al., 2000b). The full 2.5- 
200/um spectra of these three galaxies are shown in Fig. 1, and display 
the richness of the infrared range and the differences between these 
active galaxy templates (Sturm et al., 2000). The composite source 
Circinus displays broad mid-infrared features, fine structure lines of 
low excitation, H-rccombination lines and a thermal continuum that 
predominantly trace the starburst component are intermingled with 
the strong VSG and mid-infrared AGN continua and high excitation 
fine structure lines that are predominantly caused by the hard radiation 
field of an AGN. In addition, warm and cold molecular gas and warm 
atomic gas are probed by rotational lines of H2, molecular absorp- 
tion features (e.g. OH, CH, XCN, and H2O) and very low excitation 
(il3.6eV) fine structure lines (e.g. [01], [CH]), respectively. 

2.2.1. Unidentified Infrared Bands 

The mid-infrared spectra of many active galaxies are dominated by 
characteristic broad emission features, the most prominent of which arc 
located at 3.3, 6.2, 7.7, 8.6 and 11.3/xm. They are invariably observed 
together and are commonly referred to as unidentified infrared bands 
(UIBs). These features have been detected in a wide range of celestial 
sources (e.g Allamandola et al., 1995; Cesarsky et al., 1996a; Cesarsky 
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Figure 1. Composite SWS+LWS spectra of Circinus, M82 and NGC 1068 displaying 
the spectral differences between the different source types (Sturm priv. comm.). 
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et al., 1996b; Verstraete et al., 1996; Mattila et al., 1996; Lemke et al., 
1998; Mattila et al., 1999) and are commonly detected in the spectra 
of normal galaxies. UIBs may originate from C-C and C-H bending 
and stretching vibrational modes in aromatic hydrocarbons (Leger and 
Puget, 1984; Allamandola et al., 1985; Puget and Leger, 1989) and 
the favoured carriers arc polycyclic aromatic hydrocarbons (PAHs). 
The features are attributed to PAH clusters (Boulanger et al., 1998b), 
while individual PAHs can contribute to the 5-10)um continuum and 
may be responsible for the non-linear continuum rise seen in the far 
UV (Siebenmorgen and Kruegel, 1992). Alternative carriers arc small 
amorphous carbon grains that emit when exposed to moderately in- 
tense UV/visible radiation (Puget and Leger, 1989; Boulade et al., 
1996; Uchida ct al., 1998; Pagani et al., 1999). Since the exact molecular 
structure of the carriers remains under debate, the fcatTircs are referred 
to as UIBs throughout this text (for a more thorough review on the 
structure of UIBs see the chapter of Peeters et al. in this volume or 
Peeters et al. (2003)). Where strongly detected, the carriers are likely 
to provide the majority of the flux in the 3-12//m band and show near- 
invariant profiles in a variety of galaxies (Rigopoulou et al., 1999; Helou 
et al., 2000; Forster Schreiber et al., 2003). The widths of the lines 
are thought to arise from temperature broadening of the vibrational 
bands and have been fit using Lorentzian rather than Gaussian profiles 
(Boulanger et al., 1998b; Verstraete et al., 2001) but see the discus- 
sion in Peeters et al. (this volume) for further details on this issue. 
The spatial correlation between mid-infrared emission with atomic and 
molecular hydrogen suggests that UIB emission originates at the in- 
terface of HII and molecular clouds, the so-called photo-dissociation 
regions (PDRs, e.g. Verstraete et al., 1996; Cesarsky et al., 1996b; 
Pagani et al., 1999). 

Mid-infrared starburst spectra are generally dominated by these UIB 
features and show a strong continuum from very small grains. The 
relative strengths of the features are sensitive to local radiation fields, 
distribution of grains throughout the galaxy system and extinction 
(Verstraete ct al., 1996; Sturm et al., 2000; Lutz et al., 1998b; Laurent 
et al., 2000). UIB emission in starburst regions originates from UV- 
heated ionised UIB carriers. A significant contribution can also arise in 
heating of neutral UIB carriers in the extended regions by the diffuse 
ISRF i.e. by less energetic (visible) photons, which is mostly seen in 
normal galaxies and is proposed in some active galaxies (Mattila et al., 
1999; Smith, 1998; Pagani et al., 1999; Haas et al., 2002). 

In intense radiation fields (~10^ times that of the local ISRF) such 
as in metal poor environments or in the vicinity of super star clus- 
ters or AGN, UIBs are weak or absent (Contursi et al., 2000; Thuan 
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et al., 1999; Rigopoulou et al., 1999; Sturm et al., 2000; Santos-Lleo 
et al., 2001). Small grains and UIB carriers arc susceptible to high 
temperature fluctuations and photons with sufficient energy may cause 
grain sublimation or destruction. This is reflected in the UlB-to-mid- 
infrared continuum strength ratio (7|Ltm/15/xm flux density) that is 
seen to be lower (>0.5) in active regions and galaxies than in more 
quiescent regions such as in normal galaxies and the disks of some 
active galaxies [ratio~l with a large dispersion, see Fig. 1 of Vigroux 
et al. (1999) and Fig. 12 in Dale et al. (2000)]. Weak detections of UIBs 
in galaxies with AGN probably originate from the cooler material such 
as circumnuclear starburst regions or diffusely heated extended regions 
rather than from the immediate vicinity of the AGN torus. In nearby 
Seyferts, spatially resolved mid-infrared spectroscopy suggested that 
the absence or suppression of UIB emission is due to the fact that the 
dust is predominantly heated by processes related to the central AGN 
[e.g. in NGC 1068 (Le Floc'h et al., 2001), Circinus (Moorwood, 1999), 
NGC4151 (Sturm et al., 1999), Mrk 279 (Santos-Lleo et al., 2001)]. This 
was conflrmed through a comparison of ISO and high spatial resolu- 
tion ground based data of AGN and starbursts, that demonstrated the 
absence of UIBs in the nuclei of AGN-hosting galaxies (Siebenmorgen 
et al., 2004b). For unresolved sources the VSG continuum may be so 
strong that the UIB features are diluted by the over-powering contin- 
uum resulting in their non-detection (Lutz et al., 1998b). Thus at high 
redshifts, where low-metallicity- and/or AGN-dominated- systems may 
be more prevalent, the intrinsic weakness of UIB features seen in local 
systems of these types would imply that the use of UIB features as a 
tracer of similar systems at high redshift may be problematic. 

2.2.2. Fine structure lines 

Fine structure lines are tracers of nebular conditions such as excitation 
and the intrinsic ionising spectrum. By virtue of their differing ex- 
citation potentials and critical densities they provide an insight into 
the energetics and chemical composition of the regions from which 
they originate. Starburst galaxies commonly show low excitation fine 
structure lines (e.g. [Arll], [Arlll], [Fell], [Nell] and [SIII]) that have 
ionisation potentials between 13.6eV and ~50eV and mostly arise from 
HII regions that have been photoionized by their central massive stars. 
A PDR or shock origin may also contribute to the flux of some low- 
excitation lines. High excitation lines with potentials exceeding 50eV 
(e.g. [NeV,VI], [SIV],[MgV,VII,VIII] [OIV], [SiV,VI,VII]) are present 
in the spectra some low-metallicity galaxies but are more commonly 
found in AGN where high energy photons, produced by non-thermal 
processes related to the central black hole, photoionise the gas. Ex- 
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citation to these high energy levels (>50eV) is difficult to attain by 
photons emitted by stars in typical HII regions. However, the mid- 
infrared spectra of some starbursts do show the high excitation [OIV] 
and [SIV] lines. The presence of these lines is associated to collisionally 
excited, shocked or coronal gas (e.g. Lutz et al. (1998a)). Active galaxy 
spectra also show fine structure lines with excitation potentials below 
the ionisation potential of hydrogen (e.g. [OIII]52,88/xm, [OI]63,145/xm 
and [CII]158//m), that trace mostly cool atomic and molecular clouds. 
[CII]158/im and [OI]63/Ltm are important cooling lines in PDRs but also 
originate in widespread atomic gas (sec section 4.6). 

As infrared fine structure lines are fairly insensitive to uncertain- 
ties of the effective temperature of the gas, and because they are less 
affected by extinction than optical lines, they provide an unequalled 
insight into the properties of the ISM (e.g. electron temperature and 
density, ISRF, metallicity) of active galaxies and their embedded power 
source. Ratios of two lines of the same species provide information on 
the nebular conditions such as electron temperatures and densities. 
Ratios of high to low excitation states of the same element probe 
the hardness of the ISRF and ratios of high- to low-excitation lines 
of different elements have been used as tracers of AGN presence and 
thus contribute to nebular diagnostics (see section 6.1). Spectroscopic 
data from ISO was a major advance over IRAS enabling the quantita- 
tive estimation of the dominant fuelling mechanism in a given source 
(Genzel et al., 1998; Laurent et al., 2000; Sturm et al., 2000). 

2.2.3. Molecular Features 

Molecular features trace the cold dense molecular material in active 
galaxies. Several active gas-rich starbursts, Seyferts and ULIGs display 
a range of molecular lines: the rotational lines of molecular hydrogen 
(Sect. 2.2.3.2), OH, H2O, NH3 and CH are seen in absorption or emis- 
sion (see Fig. 1 of Fischer et al., 1999b, also refer to Hur et al., 1996; 
Fischer et al., 1997; Bradford et al., 1999; Fischer et al., 1999a; Fischer 
et al., 1999b; Colbert et al, 1999; Sturm et al., 2000; Rigopoulou et al., 
2002; Spoon et al., 2002; Spinoglio et al., 2005). Molecular absorption 
features appear to be more common (and prominent) in sources of high 
luminosity (e.g. ULIGs) than those of lower luminosity (Fischer et al., 
1999b; Fischer et al., 1999a; Spoon et al., 2002). Spectra displaying 
molecular absorption features show a wide variety of excitation: from 
starburst galaxies like M82 which display features indicating that most 
of the O and C molecules occupy the ground-state, to objects such as 
ULIGs where a significant population exists at higher energy levels (Fis- 
cher et al., 1999b; Fischer et al., 1999a). Additional absorption features 
at 6.85 and 7.25/xm are seen in the spectra of some active galaxies and 
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ULIGs (Spoon et al., 2002). They are attributed to CH-deformation 
modes of carbonaceous material because of their similarity to features 
seen along lines of sight towards Sgr A* (Chiar ct al., 2000). ULIG 
and Seyfert spectra also exhibit an absorption feature at 3.4//m that is 
likely to be aliphatic CH stretch absorptions of refractory carbonaceous 
material (Imanishi, 2002; Dartois et al., 2004). 

2.2.3.1. Silicate Absorption The 9. 7/nm and 18//m absorption features 
arise from the stretching and bending modes of SiO. The extinction 
suffered by a galaxy can be estimated by the depth of these lines. 
For example, the extreme 9.7/xm silicate absorption feature seen in the 
starburst-ULIG Arp 220 provides a lower limit to the optical depth of 
Ay ^45 (Charmandaris et al. (1997), but see Spoon et al. (2004) for a 
more detailed analysis of this system). Though it is possible that the 
depth of the 9.7/im may be overestimated in the feature-dense mid- 
infrared spectra of starbursts exhibiting UIB emission, (Sturm et al., 
2000; Forster Schreiber et al., 2003). The extreme extinction of the 
starburst M82 (Ay =15-60) derived by Gillett et al. (1975) was ascribed 
to this problem as ISO data could not support such a high extinction 
(Sturm et al., 2000). 

2.2.3.2. Molecular hydrogen and warm molecular gas Rotational tran- 
sitions of H2 that originate in dense molecular media were detected in 
the spectra of several active galaxies: NGC3256 (Rigopoulou ct al., 
1996); NGC891 (Valentijn and van der Werf, 1999); NGC 1068 (Lutz 
et al., 2000b); NGC 4945 (Spoon et al., 2000); NGC 6240 (Lutz et al., 
2003). Low rotational transitions originate from warm gas (T~100- 
200K) and the higher and rovibrational lines from much hotter gas 
(T^IOOOK). Rigopoulou et al. (2002) investigated a sample of starbursts 
and Seyferts observed by ISO presenting pure rotational lines from 
S(7) to S(0). Irrespective of starburst or Seyfert type, temperatures of 
~150K arc derived from the S(1)/S(0) ratio. A similar temperature 
was found for NGC 4945 (Spoon et al., 2000). While ~10% of the 
galactic gas mass in a starburst system can be accounted for by the 
warm component, on average higher fraction (18%, range 2-35%) are 
observed for Seyferts (Rigopoulou et al., 2002). The temperature of 
the molecular hydrogen implies that it arises from a combination of 
emission from fairly normal PDRs (at least for starbursts) , low velocity 
shocks and X-ray heated gas around the central AGN. The latter is 
more probable in Seyfert galaxies which accounts for the larger fraction 
of warm gas found in these systems. Shock heating appears to be the 
most likely origin of the extremely strong rotational H2 lines detected 
in the merging 'double active nucleus' system NGC 6240 (Lutz et al.. 
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2003). While, the origin of shock-heating can be largely associated to 
winds originating from massive stars and supernovae, recently Haas 
et al. (2005) have found the first observational evidence of pre-starburst 
shocks in the overlap region of the early merger the Antennae (see 
Section 4.1). Prom a re-analysis of archival ISOCAM-CYF data, they 
find that the strongest molccTilar hydrogen emission (traced by the 
H2S(3) line) is displaced from the regions of active star foramtion. This, 
together with the high line luminosity (normalised to the far-infrared 
luminosity) indicates that the bulk of excited H2 gas is shocked by the 
collision itself. 

2.2.3.3. Ice absorption and cold molecular gas An interesting discov- 
ery made by ISO was the first extragalactic detection of absorption 
features due to ices (H2O, CH4 and XCN) present in cold molecular 
components of starbursts, Seyferts and predominantly ULIGs. Sturm 
et al. (2000) first reported the detections of weak water ice absorption 
features at 3/Ltm in the spectra of local starbursts M 82 and NGC 253. 
Very strong water ice, as well as the first detections of absorptions due 
to CO and CO2 ices, were detected in the actively star-forming galaxy 
NGC 4945 (Spoon et al., 2000). The presence of ices is related to the 
presence of cold material and the radiation environment. One would 
expect cold molecular clouds in starbursts to be likely hosts while more 
intense environments, such as those predominantly fuelled by AGN, to 
be less favoured. This is corroborated by the lack of absorption features 
in the spectrum of NGC 1068 (Sturm et al., 2000) and the limits on the 
strength of the CO absorption feature set from ISOSWS spectra of 
31 AGN (Lutz et al., 2004b). However, water ice absorption has been 
detected in Seyfert galaxies, e.g. in the UIB-free spectrum of the Seyfert 
galaxy NGC 4418 (Spoon et al., 2001). The absorption dominated spec- 
trum bears strong similarities to the spectra of embedded protostars, 
and the depth of absorption features implies that its Seyfert nucleus is 
so deeply embedded that ices are shielded from the intense radiation 
field of the AGN. 

In a heterogeneous sample of 103 active galaxies with high signal- 
to-noise mid-infrared spectra, approximately 20% display absorption 
features attributed to ices (Spoon et al., 2002). While ices are weak 
or absent in the spectra of starbursts and Seyferts, in ULIGs they are 
strong, considerably stronger than the weak features found previously 
in M 82 and NGC 253. This implies that rather than the radiation envi- 
ronment it is the amount of cold material that determines their presence 
in these most luminous and massive infrared galaxies that are known 
to contain large concentrations of molecular material (e.g. Solomon 
et al. (1997); Downes and Solomon (1998)). Although plausible, more 
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sensitive (and spatially resolved) mid-infrared spectroscopy is required 
to understand these differences, and the distribution and energetic con- 
ditions of cold molecular material. Moreover, higher spectral resolution 
over a wider mid-infrared wavelength range is necessary for accurate 
decomposition of the spectra into their various components, which will 
aid identification of subtle and/or blended features (e.g. Spoon et al., 
2004). Spoon et al. (2002) proposed an evolutionary sequence for their 
ice spectra ranging from a highly obscured beginning of star formation 
where ices dominate the spectra, to a less obscured stage of star forma- 
tion where the UIBs become stronger. At any stage the star formation 
may be accompanied by AGN activity. 

2.2.3.4. OH features and (Mega- Jmasers Transitions of the OH molecule 
were also detected in the ISOSWS and ISOLWS spectra of some active 
galaxies (Fischer et al., 1999b), such as in the starbursts NGC253 
(Bradford et al., 1999) and NGC4945 (Brauher et al., 1997) permitting 
estimation of the OH column density and abundance. The depths of the 
OH features detected in the ISOLWS far-infrared spectra of galaxies 
hosting known OH (mega-)masers were found to be sufficient to provide 
all of the infrared photons required to radiatively pump the maser ac- 
tivity (e.g. in the starburst-ULIG Arp 220 (Skinner et al., 1997; Suter 
et al., 1998); AGN-ULIG Mrk 231 (Suter et al., 1998); the merging 
double systems ULIG IRAS 20100-4156 and LIG 3 Zw 35 (Kegel et al., 
1999); the starburst NGC4945 (Brauher et al., 1997)). In Arp220, the 
fact that the far-infrared continuum source is extended and that the 
maser must be located in front of it favours the interpretation that a 
starburst provides the continuum photons for its maser (Skinner et al., 
1997). The high pump rate determined for this source suggests that 
pumping mechanisms other than radiative pumping (such as collisional 
pumping) may contribute (He and Chen, 2004). 

2.3. Magnetically Aligned Dust Grains 

Polarised emission detected from active galaxies is ascribed to dichroic 
absorption by the ISM of a galaxy i.e. by elongated dust grains that 
are magnetically aligned. This has been observed in galactic sources 
where the position angle of the measured polarisation indicates the 
orientation of the projected magnetic field. Polarisation is an important 
probe of the physical conditions in active galaxies as it can provide 
constraints on the formation of the putative torus (Alexander et al., 
1999b). Such polarisation studies were performed for some ULIGs and 
revealed 3-8% of the 5-18/um flux to be polarised (Siebenmorgen and 
Efstathiou, 2001). The archetypal starburst-like ULIG Arp220 has the 
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lowest polarisation while the warm AGN-like Mrk 231 ULIG the high- 
est. The first polarisation profile of a starburst galaxy was made at 
Gfim for NGC 1808 showing a 20% increase towards the outer regions 
(Siebenmorgen et al., 2001). This, together with a complementary mea- 
surement at 170/xm (unresolved), are best explained by the large scale 
(~500pc) magnetic alignment of (non-spherical) large grains (> lOOA). 

3. Very Cold dust and diffuse radiation fields: ISM 

ISOPHOT-data clearly confirmed the presence of very cold dust (T~10- 
20K) in a variety of infrared galaxies. Often in conjunction with sup- 
plementary (sub-)mm photometry, this component was detected in 
spatially resolved maps or intensity profiles of local galaxies (e.g. Haas 
et al. (1998b)), and was isolated as a significant contributor to the 
far-infrared luminosity in SED analyses of unresolved active galaxies 
(Calzetti ct al, 2000; Klaas et al., 2001; Spinogho et al., 2002; Haas 
et al., 2003). Although cold 'cirrus' or 'disk' components in galaxies 
were indicated by IRAS colours, IRAS was insensitive to this very cold 
component due to its shorter wavelength coverage (< lOO^um). As a 
result gas-to-dust ratios of distant galaxies based solely upon ffi^S'data 
were overestimated whereas the addition of ISO measurements resulted 
in values which are closer to that of the Galactic value (~160) (Haas 
et al., 1998b) The very-cold component often accounts for a substantial 
fraction of far-infrared emission and is most commonly seen in early 
type and disk galaxies (see Sauvage et al. this volume). With a scale- 
length in excess of the stars, less than the atomic HI gas and similar 
to that of the molecular gas, it is likely that the dust is associated to 
this cold gas components (e.g. Davies et al., 1999; Popescu and Tuffs, 
2003; Popescu et al., 2004). 

Very-cold dust is also seen in both starburst and AGN-hosting active 
galaxies. NGC 253, that hosts a nuclear starburst, has been shown to 
contain extended cold dust along the major and minor axes. While the 
dust along the minor axis is attributed to outflows, the cold emission 
along the major axis is shown to have a scale length 40% greater than 
that of the stars (Radovich et al., 2001). The correlation of the 90/im 
emission profile of four Seyfert galaxies with the R-band extent but 
not with the star formation tracer Ha suggests that a large contri- 
bution of the far-infrared emission from Seyfcrts does not arise from 
the (nuclear) star-forming regions but from a cooler far-infrared disk 
component within which normal stars heat dust grains (Perez Garcia 
et al., 2000). The detection of substantial cold far-infrared emission 
(T<20 K) from the star-forming obscured overlap region of the merging 
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Antennae galaxy (see Section 4.1) indicates the presence of cold dense 
dust within which stars are forming (Haas et al., 2000a). 

SED analysis using black or grey-body decomposition suggested the 
need for more than one dust component, including a very cold dust 
component, to explain the emission above 100/Ltm of many galaxies. 
However, by allowing both temperature and dust grain emissivity to 
be free parameters, it is possible to approximate the entire far-infrared 
SED (A >50/im) using a a single modified black body. Invariably for 
active galaxies this was achieved with temperatures ~30-50K and emis- 
sivity values of 1.5-2 (e.g. Colbert et al., 1999; Siebcnmorgcn ct al., 
1999a). However, by doing so raised conflicts with other observational 
constraints (see Klaas et al. (2001) for a full description). Moreover 
by leaving the dust emissivity as a free parameter presupposes that 
intrinsic dust properties differ from source to source with no phys- 
ical basis supporting this assumption. By choosing a Galactic value 
of 7=2 for the grain emissivity (i.e. assuming that dust properties 
do not vary between sources), single black-bodies fail to explain all 
the far- infrared emission - additional cold components (T~10-20K) are 
required that are associated with the spatially resolved extended very 
cold components described above. 

The SEDs of active galaxies are warmer than those of normal galax- 
ies and the far-infrared part is generally fit by two components repre- 
senting active star formation and very cold (diffusely heated) material. 
The very cold dust can be a significant contributor in terms of galaxy 
mass and energy output. For example, two of the three components used 
to fit the SED of the LINER NGC 3079 have very cold temperatures 
of 12 and 20K, respectively (Klaas and Walker, 2002). Calzctti ct al. 
(2000) find that two dust components are required at T=40-55K and 
T~20-23K to explain the far-infrared SEDs of eight local starbursts. 
The latter comprises up to 60% of the total flux and the associated 
mass is as much as 150 times that of the former. For the most actively 
star forming galaxies or those with strongest radiation fields, the contri- 
bution to the bolometric luminosity from cold emission (as measured in 
the 122-1100//m range) was shown to be smaller in actively star forming 
(5%) than in cirrus dominated galaxies (~40%) (Dale ct al., 2001a). 
SED analysis of a wide range of active galaxies indicates that very cold 
components are common whether they are centrally fuelled by AGN or 
starbursts, suggesting that the central activity does not directly affect 
the far-infrared emission. 
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4. Cold dust and low excitation radiation fields: starburst 

regions 

The T=40-45K temperature component found in local starbursts by 
Calzetti et al. (2000) mentioned above represents the emission from 
dust surrounding sites of active star formation. For starburst galaxies 
this emission (peaking between 60 and lOO/^m) dominates the infrared 
SED and is the major contributor to the infrared luminosity. This 'cold' 
infrared emission is generally less extended than the very cold compo- 
nent described in the preceding section. In some nearby galaxies, the 
far-infrared morphology traces star forming regions well. For example, 
the strongly far-infrared-emitting dark molecular cloud coincides with 
a cluster of star forming regions in the galaxy NGC 4051 (Perez Garcia 
et al., 2000). Similarly, far- infrared emission could be resolved for the 
pair members of interacting spirals where enhanced star formation has 
been established (Xu et al., 2000; Gao et al., 2001). The correlation 
of the warmer far-infrared emission with the R-band morphology of 
local Seyferts hosting compact nuclei and circumnuclear star forma- 
tion further supports that this far-infrared component is related to 
star formation activity (Rodriguez Espinosa and Perez Garcia, 1997). 
Moreover, far-infrared emission has also been spatially associated with 
regions exhibiting signs of strong star formation (e.g. strong UIBs). 
Roussel et al. (2001) demonstrated that the integrated far- infrared flux 
and the star-formation related mid-infrared flux are strongly corre- 
lated in the disks of normal star forming galaxies indicating that the 
grains responsible for the mid- and far-infrared share a common heating 
source. 

4.1. The importance of interactions/mergers 

Interactions and mergers enhance starburst activity, both nuclear and 
extra-nuclear (see review by Struck, 1999). The final merger stage plays 
a critical role in how active galaxies form stars and may influence the 
appearance of AGN and the onset of the ULIG phase. The fraction of 
disturbed or interacting systems in infrared selected samples of galaxies 
appears to increase with luminosity with nearly 100% of all ULIGs 
displaying evidence of interaction (see Sanders and Mirabel, 1996 and 
references therein). Numerical simulations have also helped in estab- 
lishing that interactions and mergers play an important role in the 
formation and evolution of such galaxies (Mihos and Hernquist, 1996). 
The induced gravitational instabilities form bars which strip in-falling 
gas of angular momentum and enable radial inflows to the nuclear 
regions that can feed starbursts or AGN (Combes, 2001). Moreover, 
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Figure 2. This composite figure of the Antennae is from Mirabel et aL (1998) show- 
ing the ISOCAM-CYF contours overlaid on a HST V and I band image. Half of the 
mid-infrared emission arises from starburst activity that is completely obscured in 
the optical, including the brightest mid-infrared knot (knot A). The CVF spectrum 
of this knot is shown in the lower part of the image displaying strong signatures of 
massive star formation - the Nelll line and a strong rising continuum above lO^m. 
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mergers appear to be the triggering mechanism for the ultraluminous 
phase. Dynamical shocks and tidal forces are the only mechanisms 
that can cause transference and compression of sufficient quantities 
of matter to fuel luminosities in excess of 10^^ Lq. 

ISOCAMs spatial scale was sufficient to map the morphologies of 
local interacting active systems (see Mirabel and Laurent, 1999 for a re- 
view) and agree with results from ground-based high resolution images 
that (optically-faint or invisible) compact sources located in deeply 
obscured nuclear regions or interaction interfaces dominate the mid- 
infrared emission (e.g. Soifcr et al., 2000; Soifcr ct al., 2001; Charman- 
daris et al., 2002; Gallais et al., 2004; Siebenmorgen et al., 2004b). For 
some sources, diffuse emission between these 'hot-spots' can also con- 
tribute significantly to the total mid-infrared emission (e.g. Le Floc'h 
et al., 2002; Gallais ct al., 2004). This diffuse emission can be enhanced 
in merging systems that are common among active galaxies. 

By probing through the famous characteristic dust lane of the AGN- 
hosting elliptical galaxy Centaurus A (NGC 5128), imaging with ISO- 
CAM revealed a barred 'mini-spiral' located in its central 5 kpc. The 
mini-spiral was proposed to be tidal debris from a gas rich object 
that was accreted in the last gigayear (Mirabel et al., 1999; Block and 
Sauvage, 2000). Evidence of accreted companions in Cen A exists from 
previous studies of its stellar and gas kinematics (Charmandaris et al., 
2000, and references therein). /5'OPiifOT' contributed to the accretion 
scenario through the detection of the northern cold dust shell that is 
possibly the ISM remnant of a captured disk galaxy (Stickel et al., 
2004). As a result Cen A is a prime example of a giant elliptical that 
must have experienced several minor mergers. It is likely that, as the 
interactions took place, tidal forces funnelled gas into the dynamical 
centre and boosted the star formation rates in the gas rich components 
forming the dust observed in the infrared and visible. 

4.1.1. Nuclear and extra-nuclear starhursts 

Disturbed morphological signatures such as rings, tidal tails, bridges 
and arcs are clearly detected in ISOCAM observations of numerous 
merging galaxies (e.g. the Antennae (Mirabel et al., 1998) Fig. 2; UGC12915, 12914 
(Jarrett et al., 1999); NGC 985 (Appleton et al., 2002); Arp 299 (Gallais 
et al., 2004)). The mid-infrared morphologies depend upon the mass 
ratio of the progenitor galaxies and the merger stage. As well as enhanc- 
ing star formation in the nuclear regions, interactions cause a general 
increase in the widespread activity of the parent galaxies. At the early 
stages of interaction, wide scale star formation can be induced and dust 
grains arc heated by a generally warmer radiation field. For example, 
ISOCAM observations of the VV 114 galaxy show that 60% of the 
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mid-infrared flux arises from a diffuse component that is several kpc in 
size (Le Floc'h et al., 2002). Nearly forty percent of the 7 and 15/im 
emission detected in the early-stage merger (nuclei separation 5kpc) 
galaxy pair Arp 299 (NGC 3690 and IC 694) is diff'use and originates 
from the interacting disks (Gallais et al., 2004). In more advanced merg- 
ers, shocks and gravitational instabilities encourage the onset of star 
formation. The transfer of matter to nuclear regions and compression 
in interaction zones gives rise to enhanced mid-infrared emission in 
nuclear and often extra-nuclear sites. A clear example of extra-nuclear 
star formation was uncovered by ISOCAM in the intra-group medium 
of Stephan's Quintet (Xu et al., 1999). 

The best example of ISOCAMs capability in probing through highly 
cxtincted regions is imaging of the spectacular major merger in the 
Antennae (Arp 244, NGC 4038/4039). The completely obscured overlap 
region in the HST-WFPC2 image {Ay ~ 35mag Verma et al., 2003) 
was shown by ISOCAMto contain a massive star forming region (^50pc 
in radius) comprising more than 15% of the total 12-17;Um luminosity 
which outshines the two nuclei by a factor of 5 (Mirabel et al., 1998). 
The warm mid-infrared emission traces the molecular gas that fuels star 
formation (Wilson et al., 2000; Wilson et al., 2003) whereas the cold 
(T ~ 30K) and very cold dust (T 20K) are traced by far-infrared 
ISOPHOT and sub- mm observations (Haas et al., 2000a). Haas et al. 
(2000a) find that the sub-millimetre knots, in particular K2, are likely 
to be in a pre-starburst phase and the simultaneous presence of power- 
ful off-nuclear starbursts &; pre-starbursts may be a general feature of 
colliding galaxies. The authors postulate that once star formation has 
commenced in these clouds, the Antennae may evolve from a luminous 
IR galaxy into an ultraluminous one. 

An interesting feature for the Antennae (also seen for Centaurus 
A) is that the most luminous 15/Ltm knots are not coincident with 
the darkest, most prominent dust absorption lanes seen in the HST 
image. Presumably the darkest optical dust lanes trace the coldest 
emission (i.e. that probed by the far-infrared) rather than the warm 
dust from which 15/xm emission arises. This is plausibly the cause of 
the displacement. Alternatively projection effects may be responsible 
for this offset (Mirabel et al., 1998). 

Extranuclear overlap starbursts emitting in the mid-infrared are 
not unique to the Antennae. In a sample of 8 spiral-spiral pairs with 
overlapping disks, Xu et al. (2000) found extranuclear, overlap re- 
gion starbursts in 5 galaxies classified as advanced mergers or having 
severely disturbed morphologies (but not in less disturbed systems). 
However, the starbursts in the bridges or tidal tails are generally fainter 
than those in the nuclear regions. 
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4.1.2. Collisional Rings 

Collisional ring galaxies arc produced when a compact galaxy passes 
through the centre of a larger disk galaxy giving rise to a radial density 
wave that creates a massive star forming ring (see review by Appleton 
and Struck-Marcell, 1996). If the collision is off centre then the ring 
morphology is asymmetric e.g. as in VII Zw 466 (Appleton et al., 1999). 
These systems are excellent local examples of induced star formation 
in colliding galaxies. The effects of collisions on the ISM can be in- 
vestigated in detail, and may be relevant for higher redshift galaxies 
including ULIGs (Appleton et al., 2000). Azimuthal variations in over- 
density in the rings are often noted as well as extranuclear star forming 
knots that show strong UIBs (Charmandaris et al., 1999; Appleton 
et al., 1999; Charmandaris et al., 2001b; Appleton et al., 2002). When 
detected, star formation in the ring contributes at a level of Z 10% of 
the bolometric output (Appleton et al., 1999; Appleton et al., 2000). 
The mid-infrared emission is clearly enhanced in the nucleus and the ex- 
panding ring even though diffuse emission from the intervening medium 
was also observed (Charmandaris et al., 1999; Charmandaris et al., 
2001b). The diffuse components coincide with regions where only weak 
Ho; had been previously measured. The mid-infrared colours measured 
in collisional ring galaxies are similar to those of late-type galaxies, 
indicative of more modest star formation than is found for extranuclear 
starbursts (e.g. in the nucleus inner ring and spokes of Cartwheel, the 
entirety of ArplO and the ring of Arp 118 (Charmandaris ct al., 2001b)). 
The morphological agreement between the mid-infrared and Ha or ra- 
dio emission is good although differences have been noted for VII Zw 
466 (Appleton et al., 1999). A deviation from this picture is seen for 
the extranuclear starburst knot in the archetypal collisional ring galaxy, 
the Cartwheel, that coincides with the strongest Ha and radio knot. 
The regions has an unusually red mid-infrared colour (15/7/Ltm~5.2) 
which is the highest among the extragalactic regions observed with 
ISO CAM and it is even more extreme than the powerful starburst 
region (Knot A) revealed in the Antennae (Charmandaris et al., 1999). 
Interestingly this region also hosts the remarkable hyperluminous X- 
ray source (Lo.s-iofceV ~ 10^^ ergs/s Gao et al., 2003) which dominates 
the X-ray emission from the Cartwheel ring. 

4.2. Hot star population, Starburst Evolution & 
Metallicity 

Starburst galaxies, defined as having star formation rates that cannot 
be sustained over a Hubble time, arc likely the birthplace of a large 
fraction of massive stars. Locally, four starburst galaxies can account 



Page 20 



Obscured "Activity": AGN, Quasars, Starbursts and ULIGs 



21 



for 25% of the local massive star population within lOMpc (Heckman 
et al., 1998) and arc important as the providers of the Lyman con- 
tinuum photons that ionise hydrogen. 75*0 spectroscopy enabled the 
stellar populations of starbursts to be probed using fine structure lines 
that are sensitive to the nebular conditions from which they originate 
(see Forster Schreiber et al. (2001) for a detailed analysis of the stellar 
populations and radiation environment of the local starburst M82). 
The ratio of [NeIII]15.6/xm (Ep=41eV) and [NeII]12.8/xm (Ep=22eV) 
is sensitive to the hardness of the stellar energy distributions of the OB 
stars which excite them. In a survey of 27 starburst galaxies, Thornley 
et al. (2000) confirmed earlier ground-based results that the excitation 
measured by the [NeIII]15.6/im/[NcII]12.8//m ratio was lower in star- 
burst galaxies than measured for Galactic HIT regions, indicating a 
deficiency of massive stars in the starbursts. Thornley ct al. suggest a 
lower upper mass cut-off to the initial mass function (IMF) or ageing 
of the massive star population to be responsible for the low excitation. 
Through extensive quantitative photo-ionisation modelling, Thornley 
et al. found the ratios to be consistent with the formation of mas- 
sive stars (5O-IOOM0). This, together with the known presence of very 
massive stars in local starburst regions, led Thornley et al. to prefer 
an ageing scenario to explain the low excitation. This explanation was 
further supported by modelling presented in Giveon et al. (2002). The 
modelling of Thornley et al. indicates a short starburst timescale of only 
a few O star lifetimes (lO^-lO'^ yr) and suggests that strong negative 
feedback regulates the star-forming activity. Disruption due to stellar 
winds and SNe may cease the starbursts before the depletion of gas and 
suggests that periodic starburst events are likely (e.g. in M82 Forster 
Schreiber et al. (2001)). 

Verma et al. (2003) confirmed that for a given nebular abundance, 
local starburst nuclei are of lower excitation than Galactic and Mag- 
ellanic cloud HII regions. This was shown for the primary products 
of nucleosynthesis argon, sulphur and neon confirming the Thornley 
et al. (2000) result was not restricted to Ne only. Additionally, the 
observed excitation was confirmed to be metallicity dependent, with 
the lowest metallicity compact dwarfs showing the highest excitations, 
at a similar level to local HII regions. Ageing, a modified IMF and 
metallicity are all possible contributors to the metallicity-excitation 
relation. In the same study the abundance of sulphur was observed 
to be lower than expected for a primary product of nucleosynthesis. 
Wolf-Rayet and BCD galaxies did not show such an under-abundance 
(Verma et al., 2003; Nollenberg et al., 2002). This underabundance and 
weakness in the sulphur line strength implies that (a) the fine structure 
lines of neon are favoured over the those of sulphur as a star formation 
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tracer in future spectroscopic surveys of galaxies and (b) that sulphur 
is an unsuitable tracer of metallicity. 

4.3. Wolf-Rayet Starbursts 

The Wolf-Rayet phase is a signature of massive star formation as the 
progenitors are believed to be massive young O stars (3 ^ tsb $ 8Myr, 
M ^ 2OM0). WoIf-Rayct features are found in a wide range of galax- 
ies such as BCD's, HII galaxies, AGN, LINERs and starbursts (see 
Schaerer et al., 1999 for a compilation). Recently, Li'pari et al. (2003) 
identified possible Wolf-Rayet features in the spectra of some PG quasars 
The identifying Wolf-Rayet signatures are in the visible and thus are 
only representative of the unobscurcd components. Nevertheless a cor- 
relation between the star formation indicator UIB7.7/xm/15/xm con- 
tinuum with Ha in Wolf-Rayet galaxies suggests that the infrared 
and visible trace, at least partly, the same components (e.g. in Haro 
3 Metcalfe et al., 1996). The short-lived Wolf-Rayet phase is com- 
monly detected in low metallicity systems with a strong ISRF such 
as blue compact dwarfs (BCDs, see Sec. 4.4). ISO observations of 
known Wolf-Rayet galaxies reveal differences in their infrared proper- 
ties compared to non- Wolf-Rayet galaxies. In comparison to starbursts, 
Wolf-Rayet galaxies are of higher excitation and lower abundance con- 
trary to expectation of a more important role of Wolf-Rayet stars at 
higher metallicities from stellar evolution models (Verma et al., 2003). 
The Wolf-Rayet galaxy NGC 5253 has a weaker excitation determined 
from mid-infrared ISO data than expected by stellar models (Crowther 
et al., 1999) implying that the UV spectrum may not be as hard as 
commonly thought. However, Crowther et al. (1999) show that the 
effect is due to a young compact starburst surrounded by an older 
starburst both of which lie within the SWS aperture. As a result, 
measured line strengths are diluted and the apparent effective tem- 
perature of the emitting region is lowered thus reducing the measured 
excitation. This effect, that is fully consistent with an ageing/decaying 
starburst scenario suggested by Thornley et al. (2000), can explain the 
low excitation of starbursts described in Section 4.2 (Thornley et al., 
2000; Verma et al., 2003). 

4.4. Blue Compact Dwarfs 

A number of low metallicity (Z<<Zq) blue compact (<5kpc) dwarf 
galaxies were also observed by ISO. These include I Zw 36 (Mochizuki 
and Onaka, 2001; Nollenberg et al., 2002); Haro 11 (Bergvall et al., 
2000), Haro 3 (Metcalfe et al., 1996), NGC 5253 (Crowther et al., 1999), 
H Zw 40 (Schaerer and Stasihska, 1999; Madden, 2000), NGC 1569 
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(Madden, 2000; Haas et al., 2002; Lu et al., 2003; Galliano et al., 2003), 
NGC 1140 (Madden, 2000)). Super-star clusters (SSCs) have been iden- 
tified within some of them which are likely to be sites of massive star 
production that generate the bulk of the ionising flux (Gorjian et al., 
2001; Vacca et al., 2002; Vanzi and Sauvage, 2004)). These clusters can 
also account for a substantial fraction of the emerging mid-infrared 
emission (e.g. Vacca et al. (2002)). The high obscuration, hot infrared 
emission and thermal radio spectra (measured for some) imply these 
are embedded young clusters. The properties of BCDs are similar to low 
metallicity regions in the SMC (Contursi et al., 2000). Some BCDs also 
exhibit Wolf-Rayet features in their nuclear regions indicative of young 
massive stars producing a significant amount of ionising UV photons 
(see Section 4.3). ISO observations (see Madden, 2000 for a review) re- 
veal varying characteristics between BCDs although fine structure lines 
common to starbursts are present e.g. [OIV], [SIII,IV] and [NeII,III] 
with ratios that indicate high excitation (e.g. Crowther et al., 1999). 
The low metallicity and dust mass of BCDs in comparison to normal 
galaxies permits photoionisation on galaxy wide scales (Madden, 2000). 
UIBs are weak or absent due to destruction of their carriers by the 
intense radiation field and the mid-infrared is dominated by a strong 
VSG continuum. [CII]/(CO) emission (and [CIIJ/FIR to a lesser degree) 
can be greatly enhanced (xlO) in comparison to normal metallicity 
starburst galaxies (Bergvall et al., 2000; Madden, 2000). Along with 
the hot dust detected in BCDs (e.g. Metcalfe et al., 1996), very cold 
dust components were surprisingly observed in some BCD members 
of the Virgo cluster (Popescu et al., 2002) and in NGC 1569 (GalUano 
et al., 2003). In the latter, a very cold component (T=5-7K) contributes 
between 40-70% of the total dust mass of the galaxy. Molecular gas 
(CO) and infrared emission indicate that, although relatively metal 
poor, these dwarfs have high dust and gas densities for their optical 
sizes (Mochizuki and Onaka, 2001). 

An exceptional BCD is SBS 0335-052, one of the most metal poor 
galaxies currently known (I/AIZq). This dwarf harbours obscured, 
central SSCs (Dale et al., 2001b; Hunt et al., 2001) and shows a young 
thermal radio spectrum (Takeuchi et al., 2003; Hirashita and Hunt, 
2004). SBS 0335-252 is unusually bright in the mid-infrared. No UIBs 
have been detected and the [SIV] and [Nell] lines are weak in compar- 
ison to the strong continuum (Thuan et al., 1999). Estimates of the 
optical depth vary between authors from optically thin (Dale et al., 
2001b) to Ay ~10-30 (Thuan et al., 1999; Plante and Sauvage, 2002). 
The current optical, radio and infrared measurements make it difficult 
to present a consistent scenario for this intriguing galaxy. Further obser- 
vations are needed to elucidate how such as small galaxy can produce 
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such strong infrared emission without polluting its ISM with metals 
(Thuan et al., 1999). 

4.5. Star Formation Rates 

The non-linear correlation (with large dispersion) between the far- 
infrared and Ha emission in normal galaxies (Saiivagc and Thuan, 
1992) suggests the use of the far-infrared as a direct SFR estimator 
has limitations (Kennicutt, 1998). At least for normal galaxies, both a 
contribution from the nuclear regions (Roussel et al., 2001) as well as 
from cold dust that is unrelated to the starburst activity may be respon- 
sible for the non-linear behaviour. At the highest infrared luminosities, 
such as in a sample of ULIGs, this very cold component comprises 
only ^10% of the total far- infrared flux (Klaas et al., 2001) . However, 
for less luminous sources that often have more dominant very cold 
components (e.g the starbursts in Calzetti et al., 2000), the contribution 
to the far-infrared from heating from non-ionising, older stars should 
be considered when SFRs are determined (Lonsdale Persson and Helou, 
1987). A more basic problem is that SFR estimates based on the far- 
infrared flux are prone to errors arising simply from poorly determined 
far-infrared luminosities due to inadequate sampling or constraints on 
the far-infrared SEDs. This affects often the faintest most luminous 
sources detected in IRAS and ISO flux limited surveys where for some 
sources only one or two far-infrared data points are available. 

The mid-infrared also provides an alternative measure of the SFR. 
Based on ISO data for a large sample of non-AGN-hosting galaxies. 
Dale et al. (2001a) highlighted that, rather than the commonly used 
far-infrared flux (42-122/im), the mid- infrared SED between 20-42;um is 
most sensitive to the level of star formation activity. Additionally, the 
correlation of the mid-infrared flux 7^m/15;um flux ratio with optical 
tracers of star formation such as Ha in a range of galaxies (e.g. in 
normal galaxies (Sauvage et al., 1996; Roussel et al., 2001), Wolf-Rayet 
galaxies (Metcalfe et al., 1996), mixed pair local interacting galaxies 
(Domingue et al., 2003; Domingue et al., 2005), collisional ring galaxies 
(Charmandaris et al., 2001b), other interacting galaxies (Lc Floc'h 
ct al., 2002)) and the mid-infrared H-recombination lines, that are 
correlated to the Lyman continuum, shows that the mid-infrared flux 
can be used as a SFR indicator (Vigroux et al., 2001; Forster Schreiber 
et al., 2004). Forster Schreiber et al. (2004) provide an empirical cal- 
ibration of the UIB flux between 5-8.5/iim range as a SFR estimator 
which is valid for starbursts up to the rate of M82, close to which 
the ISRF causes sublimation of the UIB carriers. The VSG continuum 
represented by the monochromatic 15/xm flux, that has been shown 
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to correlate well with the [ArII]6.99/xm line (Forster Schreiber et al., 
2003), can withstand harder radiation fields providing a SFR estimator 
for higher redshift, more luminous sources. A potential scatter in this 
correlation can be introduced by an AGN contribution to the VSG 
continuum; an event that may become more common in galaxies of high 
luminosity (Sect. 6.5). In this case the mid-infrared may not be a direct 
indicator of the SFR unless the AGN contribution can be accurately 
subtracted (Laurent et al., 2000). 

4.6. FIR SPECTROSCOPY & THE [CII] DEFICIT 

The far-infrared spectrum of the starburst galaxy M82 (Colbert et al., 
1999) shows several fine structure lines including those that are com- 
monly detected in active galaxies: [OI]63,145 fiui; [OIII]52,88/um; [NII]122^m; 
[NIII]57/im and [CII]158/um. These features are superimposed on a 
strong thermal continuum. [CII]158/im, along with [OI]63,145;um and 
[NII]122/im, are important cooling lines of star forming regions and 
have been shown to correlate with the dust continuum over 4 orders of 
magnitude in the Galaxy (Baluteau et al., 2003). For many galaxies, 
especially those with low metallicities, [CII]158jLtm is the strongest far- 
infrared line, corresponding to 0.1-1% of the far-infrared flux (Malhotra 
et al., 1997; Madden, 2000; Bergvah et al., 2000). The spatial correlation 
of the [CII]158/Lim emission with obscured regions and the fact that it 
scales linearly with UIB and CO relative strength suggests it originates 
in PDRs (linger et al., 2000; Malhotra et al., 1997; Malhotra et al., 
2001; Masegosa et al., 2001). Fitting of the LWS spectrum of M82 
confirms this view with 75% of the line produced in PDRs and the 
remainder from HII regions (Colbert et al., 1999) while, for a sample 
of nearby galaxies, Negishi et al. (2001) find half to originate in PDRs 
and half in low density ionised gas. 

The spectra of active galaxies range from those exhibiting strong 
ionic fine structure emission lines (such as M82) to molecular absorption 
line dominated systems (such as Arp 220) (Fischer et al., 1999b; Fis- 
cher, 2000). The fine structure lines that are predominant in starburst 
galaxies are absent or weak in the spectra of higher luminosity (Malho- 
tra et al., 1997; Fischer et al., 1997; Fischer et al., 1999b; Fischer et al., 
1999a). In particular, the [CII]158/nm/FIR continuum ratio decreases 
with increasing luminosity (Malhotra et al., 1997) and in ULIGs is only 
10% of that measured for starbursts (Luhman et al., 1998; Fischer, 
2000; Luhman et al., 2003) and is even lower in LINERs (Sanei and 
Satyapal, 2002). As an important cooling line of warm atomic gas, 
[CII] was expected to be high for strongly star-forming galaxies with a 
high UV radiation field. Contrary to intuition, the absence of the fine 
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structure lines and the strong molecular absorption detected (e.g. OH, 
H2O, CH, NH3 and [01]) imply a weaker radiation field in ULIGs in 
comparison to nearby starbursts and AGN that is probably caused by 
shielding (Fischer, 2000). 

For sources deficient in [CII]158, the [OI]63,145 /xm lines appear to 
become the preferred means of cooling star forming gas (Malhotra et al., 
1999; Brauher and Lord, 2000). In a systematic survey of LWS spectra 
of 36 local galaxies, Negishi et al. (2001) found a strong decrease in the 
the [CII]158/im/FIR and [NII]122(Ltm/FIR cooling line ratios with bluer 
far- infrared colours. However no such trend was seen for [OI]63/FIR 
implying that the oxygen lines do not suffer from the [CII] deficiency. 

A number of explanations have been proposed for the [CII] deficit. 
Tracing [CII] towards the Sgr B2 complex, which has an intense radia- 
tion field and gas density similar to active galaxies, raises the possibility 
that self-absorption may contribute to the observed deficiency (Cox 
et al., 1999). The [CII]158/Ltm line is seen both in emission (origi- 
nating from Sgr B2) and absorption (Vastel et al., 2002) consistent 
with the presence of intervening (cold) clouds along the line of sight. 
This indicates the [CII]158//m can be optically thick. However, the 
weakness of the [OI]145/Ltm line observed in Arp 220 and Mrk 231 
and the large extinctions implied (Ay ^400, that exceed most other 
extinction measurements but see Haas et al. (2001) for Arp 220) pose 
problems for solely a self-absorption explanation of the [CII] /FIR de- 
ficiency (Luhman et al., 2003). Lower efficiency of photoelectric gas 
heating for high radiation fields, such as those in the PDRs of active 
galaxies, has been suggested (Malhotra et al., 1997) but an accurate 
analysis of this scenario remains to be performed. Additionally, the 
radiation field and gas density may play a role. Ageing bursts were 
discounted because of the detection of the high excitation lines [Nelll] 
and [SIII] and because older starbursts are unlikely to fuel such high 
luminosities (Fischer et al., 1999a). Increased collisional de-cxcitation, 
rather than the decrease in the photoelectric heating efficiency, was 
supported by Negishi et al. (2001). Luhman et al. (1998) proposed high 
UV flux incident on moderate density PDRs to explain the deficiency. 
However, in a revision of this work, (Luhman et al., 2003) preferred a 
non-PDR contribution to the far-infrared to explain the deficiency in 
the [CII] /FIR ratio. The authors suggest dust-bounded photoionisation 
regions as a possible origin since such regions contribute to the far- 
infrared (i.e. increased UV-absorption by dust) but do not generate 
significant amounts of [CII] . This interpretation is corroborated by the 
increasing fraction of dust-absorbed photons with ionisation parameter 
(Voit, 1992). This scenario is also preferred in the radiative-transfer 
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models of Arp 220, although the effects of extinction cannot be ruled 
out (Gonzalez- Alfonso et al., 2004) . 

The [CII] deficiency affects possible future searches of [CII]158/im 
line emission from high redshift sources that have similar rest-frame far- 
infrared colours to local ULIGs. Moreover, if the stronger cooling line 
in warmer infrared galaxies is [OI]63/im, then the recent non-dctcctions 
in four z=0. 6-1.4 ULIG-quasars (Dale et al., 2004) also indicate a pes- 
simistic view on either of these lines as high-redshift tracers of active 
galaxies. The [OI]63/xm line may also be self-absorbed. 

4.7. Gas-to-dust ratio and starburst evolution 

The gas-to-dust ratio is commonly used as an indicator of the fraction of 
gas that is yet to be turned into stars. Dust mass estimates determined 
from ISO SED analyses are used for this determination. The gas-to- 
dust ratios determined for starbursts, ULIGs, and quasars, are similar 
or only slightly larger (~ 150-400) than those determined for the Milky 
Way (~150) (Calzetti et al., 2000; Haas et al., 2003). The highest 
gas-to-dust ratios are seen in BCD galaxies indicating young stellar 
evolution (Dale et al., 2001b). Low gas-to-dust ratios and low dust 
temperatures indicate a post-starburst system (e.g. the radio galaxy 
MG1019+0535 (Manning and Spinrad, 2001) and HyLIG F15307+3253 
(Verma et al., 2002)). 

5. Warm dust and hard radiation fields: nuclear AGN 

5.1. AGN IN MERGERS 

ISO's low spatial scale could not discern emission from the torus and 
circumnuclear starburst regions even in the most nearby active galaxies. 
As AGN emission is unresolved in /S'OC^Af images, in this section we 
discuss the environment of the AGN, i.e. the heating of the ISM of the 
host galaxy. The collisional ring galaxies NGC985 and Arp 118 both 
host Seyfert nuclei. In NGC 985, near-infrared images resolve the bulge 
of a small, second galaxy (whose accretion likely led to the formation 
of the ring) near the active bright nucleus (Appleton et al.. 2002). The 
nuclei in both of these galaxies show a hot mid-infrared continuum 
that can wholly be attributed to dust heating by active nuclei. UIB 
emission is suppressed around the nuclei, suggesting they are absent in 
the vicinity of the AGN. In the case of Arp 118, the Seyfert nucleus 
accounts for 75% of the mid-infrared emission from the galaxy. Two 
other LIGs, W 114 and Arp 299, show strong and hot mid-infrared 
continuum emission, i.e. that hot dust is present (a 230-300K blackbody 
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describes the VSG hot continuum), indicating that embedded AGN 
may reside within these galaxies (Le Floc'h ct al., 2002; Gallais et al., 
2004). Seyfert nuclei seem to have high mid-infrared luminosities, as 
well as contributing significantly to the total mid-infrared emission of 
the system. In three merging pairs of ULIGs within which one member 
contains an AGN, the AGN hosting galaxy is invariably more infrared 
luminous than the other (Charmandaris et al., 2002). Mirabel and Lau- 
rent (1999) note that the most luminous sources are those that exhibit 
both AGN and starburst activity, reflecting the well known trend of 
increasing AGN importance with infrared luminosity (see Sec. 6). 

5.2. Warm SEDs 

ISO observations of the 12/xm IRAS galaxy sample (Rush et al., 1993), 
that consists predominantly of Seyfert 1 and 2s but includes starbursts 
and normal galaxies, displays a range of SED shapes that can be ex- 
pected for these source types (see Figure 2 in Spinoglio et al. (2002)). 
These average SEDs show clear differences between normal galax;ies, 
starbursts and the Seyferts in the optical to mid-infrared range, but 
they appear indistinguishable in the far-infrared. 

By approximating the Seyfert SEDs using black-body functions, au- 
thors such as Rodriguez Espinosa et al. (1996) or Anton et al. (2002) 
found that two components best fit their SEDs where the warmer com- 
ponent was associated to the active nucleus and the colder component 
with star forming or extended regions. Similar associations were made 
for the two (or three) components implied by a novel inverse Bayesian 
method to analyse the SEDs of 10 Seyfert galaxies (Perez Garcia et al., 
1998). Prieto et al. (2001) showed that the fiuxes of high excitation 
coronal lines detected in the ISOSWS spectra of Seyfert galaxies, which 
are ascribed to excitation by the intense radiation field of the active 
nucleus, were found to be correlated to this mid-infrared component but 
not to the far-infrared component. This implies that the mid-infrared 
component arises from dust heated by processes related to the active 
nucleus. 

5.3. lONISATION IN AGN AND SeYFERTS: THE BiG BLUE BUMP 

A thin accretion disk that fuels a black hole displays a generic signature 
(quasi-thermal photons at T~10^) in the extreme UV (EUV) commonly 
called the big blue bump (BBB). The bump is likely to extend from 
the UV towards the soft X-ray bump but, due to galactic and intrinsic 
absorption between the Lyman edge and a few hundred eV, the peak of 
the bump is not directly detectable. A novel method to reconstruct the 
BBB is to use the high ionisation potential infrared fine structure lines 
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(coronal lines, 50-300eV) that originate in the narrow-line region of 
Seyferts and quasars and directly probe the EUV SED. Reconstruction 
of the BBB was attempted using photo-ionisation modelling of the 
Seyferts NGC4151(Sey 1.5), NGC 1068 (Sey 2) and Circinus (Sey 2). 
The reconstruction is not only dependent upon the ionising continuum 
but also on the geometry of the NLR and the ionisation parameter. It 
was possible to constrain the presence of this bump for these galaxies 
by complementing the ISO data with UV, optical, and near-infrared 
lines. All three investigated sources (Moorwood et al., 1996; Alexan- 
der et al., 1999c; Alexander et al., 2000) require the presence of a 
thermal blue bump, a single power law continuum fails to reproduce 
the observed narrow line ratios. Circinus provides the clearest example 
of an EUV bump that peaks at 70eV and contains half of the AGN 
luminosity. The predicted SEDs of NGC 4151 and NGC 1068 display 
clear troughs falling sharply beyond the Lyman limit and rising sharply 
above lOOeV. The preferred explanation is that it is due to absorption 
by neutral hydrogen of column density Nh ~ 5 x lO^^cm"^. In the 
case of NGC 4151 this absorbing gas was identified through an analysis 
of UV absorption lines (Kriss et al., 1995). This scenario allows for 
the presence of (obscured) BBB in all of the sources investigated (sec 
also Prieto et al. (2001)). However, Binette et al. (1997) find that 
photoionization by a simple power law suffices to explain Circinus, if 
an optically thin (density-bounded) gas component is included. In a 
detailed metallicity analysis, Oliva et al. (1999) find that either ap- 
proach can equally explain the observations, and note that the ionising 
continuum is mostly dependent upon the gas density distribution. More 
complex geometry and distribution of gas (e.g. 3-D multi-cloud models) 
and emission mechanisms (shock heating and photoionisation) (Contini 
et al., 1998; Contini et al., 2002; Martins et al., 2003), than an assumed 
single cloud with filling factor <1 (Alexander et al., 1999c), are more 
appropriate assumptions to model the narrow-line emission properties 
surrounding AGN (Martins et al., 2003). 

While such detailed modelling was not possible for less well ob- 
served sources, a comparison of spectroscopic tracers, such as electron 
density and excitation, to photoionisation models enabled at least a 
consistency check to be performed between power-law or black-body 
continua. However, neither Sturm et al. (2002) nor Prieto et al. (2001) 
were able to distinguish the two for their sources with the range of lines 
they used. The metallicity, geometry and density of obscuring clouds 
and extinction all play a role in these analyses. 
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5.4. The dusty torus and variability 

If near and mid-infrared emission from the inner region of AGN is 
thermal re-radiation by dust surrounding the central engine then, over 
sufficiently long timcscales, variability seen in the optical/UV should 
be detectable in the infrared. This would provide convincing proof 
and constraints on the presence, temperature and distribution of dust 
grains from the variable source. However ISO monitoring of variable 
AGN, that were contemporaneous to other multiwavelength campaigns, 
failed to show any convincing thermal variability (e.g. Mrk 279 Santos- 
Lleo et al., 2001; BL Lac PKS 2155-304 Bertone et al., 2000; BL Lac 
2007+777 Peng et al., 2000; OVV quasar 3C 446 Leech et al., 2002). 
Rather, the variability where detected, was consistent with non-thermal 
components such as jets (Peng et al., 2000; Leech et al., 2002). 

These findings are consistent with variability traced by polarisation 
measurements of the OW quasar 3C 279 at 170/im (Klaas et al., 1999) 
showing that the polarisation detected was aligned to the radio axis and 
thus consistent with synchrotron emission. The rapid variability noted 
in two measurements separated by '^lyr confirms that the far-infrared 
emission is compact and originates from the core. 

5.5. Unification 

Under the standard unification scheme (Antonucci, 1993), the orienta- 
tion relative to the line of sight of a dusty putative torus surrounding 
the accretion disk of the central black hole, may connect various types 
of active galaxies: radio galaxies to quasars and Seyfert Is to Seyfert 
2s. A check of these schemes is offered through their infrared emission: 
the dust torus intercepts optical-UV photons released by the AGN and 
re-emits them in the infrared. At A>30;um the infrared emission is 
largely optically thin and hence isotropic. As a result the infrared SEDs 
should be similar for all members of a unified class, irrespective of the 
orientation of the obscuring torus to the line of sight. In this section we 
discuss how ISO results have contributed to understanding this issue. 

5.5.1. Radio galaxies and quasars 

Orr and Browne (1982) and Barthcl (1989) proposed a scheme that 
unifies quasars and radio galaxies as intrinsically similar sources, with 
observed differences ascribed to the viewing-angle with respect to the 
black hole. A quasar with a flat radio spectrum is observed with a 
viewing- or aspect-angle of 0° to the line of sight (pole-on) such that a 
beamed jet with superluminal motion causes the bright radio and far- 
infrared appearance (e.g. FR 0234-1-28). At intermediate angles between 
the pole- and edge-on orientations, a quasar with a steep radio spectrum 
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is seen. The shallower the aspect angle, the more obvious the thermal 
bump of the torus becomes. At large orientations (90° or perpendic- 
ular) to the line-of-sight, we observe the radio galaxy with the torus 
edge-on (e.g. CygA, Fornax A, Centaurus A) whereby the nucleus that 
contains an "optical quasar" is obscured by the dusty torus. Powerful 
Fanaroff-Riley (FR) II radio galaxies are " edge-on" quasars viewed at 
high inclination (Haas et al., 1998a). Overwhelmingly for samples of 
radio galaxies and radio quasars their SEDs appear indistinguishable 
in the far-infrared (Haas et al., 1998a; van Bemmel et al., 2000; Meisen- 
heimcr ct al., 2001; Andrcani ct al., 2002; Haas et al., 2004) confirming 
the expectation from unified schemes i.e. that the dust is isotropic and 
optically thin. This is contrary to indications from IRAS data that 
implied radio galaxies have lower mid-far infrared emission than radio 
quasars. The lower rcdshift range probed by IRAS may account for this 
(Meisenheimer et al., 2001). 

5.5.2. Narrow line radio galaxies; broad line radio galaxies and 

quasars 

The 3C sample of quasars selected at 178MHz contains two types 
of galaxies: (1) those exhibiting broad lines: steep spectrum quasars 
and broad line radio galaxies; and (2) those exhibiting narrow lines: 
FRII narrow line radio galaxies. Under the unification scenario, sources 
of identical isotropic lobe power should have the same far-infrared 
isotropic power. Haas et al. (2004) confirm earlier results of Meisen- 
heimer et al. (2001) for all SCR sources present within the ISO Data 
Archive^, providing clear evidence in favour of geometric unification. 
The considerable dispersion in lobe power can be attributed to environ- 
ment and evolution. Siebenmorgen et al. (2004a) confirm earlier work 
of Preudling et al. (2003) and find that the SEDs of narrow-line radio 
galaxies to be colder {Xpeak ~ lOO/xm) than the broad-line counterparts 
{Xpeak ~ 40/xm). This result is fully consistent with unification schemes 
where, in broad-line radio galaxies, the BLR is exposed and so is hot 
dust that causes the SED shape to appear warmer. 

5.5.3. FR I and FR II galaxies 

These extended radio sources are separated by the fact that FRII 
galaxies show powerful edge-brightened extended double lobes whereas 
FRIs have diffuse lobes that are brightest close to the nucleus. More- 
over, FRIs arc of lower luminosity than FRIIs. In a recent study of 
3CR sources, Miiller et al. (2004) found that the low radio power FR I 
galaxies have low mid- and far-infrared luminosities. In contrast, the FR 

^ http://www.iso.vilspa.esa.es 
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II sources are infrared luminous. The high far-infrared-to-mid-infrared 

ratio and the temperature of the far-infrared peak which is comparable 
to that of quiescent radio-quiet elliptical galaxies, implies that heating 
by the ISRF of the host galaxy is probable. This together with the 
finding that the dust masses of FR Is and FR lis are similar but with 
low gas fractions in the nuclear regions, implies that the black holes in 
FRI galaxies are fed at a lower rate than FRII galaxies (Miiller et al., 
2004). 

5.5.4. Flat spectrum and Steep spectrum Quasars 

Flat spectrum quasars show a strong synchrotron component (e.g. FR 
0234+28) , smoothly concatenating the mm flux to that observed in the 
mid-infrared (Haas et al., 1998a). A thermal far- infrared component 
may equal the synchrotron in strength (e.g. 3C 279), but in most cases 
only upper limits to a thermal contribution can be derived due to 
the dominating synchrotron emission. This prevalence of far-infrared 
synchrotron emission is consistent with Doppler boosted synchrotron 
emission of a relativistic jet seen almost end-on and does not contradict 
this unification (Haas et al., 1998a). 

5.5.5. Radio-loud and radio- quiet quasars 

A complete range of radio-loud, radio-intermediate and radio-quiet 
quasars was analysed by PoUetta et al. (2000). Except for flat radio 
sources, the mid- and far-infrared emission is of thermal origin. The far- 
infrared emission is predominantly explained by heating from the AGN 
itself and, while starburst components are present, they contribute to 
the infrared emission at different levels but always less than the AGN 
(<27%). 

5.5.6. Seyfert Is and Seyfert 2s 

Within the unified scheme, Seyfert Is and 2s differ only in viewing angle 
of the torus obscuring the central black hole. In Seyfert Is the broad 
line region is exposed whereas in Seyfert 2s it is hidden from view due 
to obscuration by the torus. The central regions where dust is heated to 
high temperatures gives rise to a mid-infrared peak in the continuum 
emission. Aspect-oriented unification is consistent with the result that 
the mid-infrared dust temperatures of Seyfert Is are higher than in 
Seyfert 2s, where the torus is partially optically thick up to the mid- 
infrared (Perez Garcia and Rodriguez Espinosa, 1998). Spectroscopic 
tracers such as high excitation lines and broad lines are indicative of 
the presence and orientation of the putative torus. In a large sample 
of Seyferts, Clavel et al. (2000) showed that UIB equivalent widths are 
larger in Seyfert 2s than in Seyfert Is. An interpretation of this result 
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is that UIB features originate from extended regions and hence emit 
isotropicahy, while the AGN continuum is absorbed by the torus and 
must emit anisotropically (Perez Garcia and Rodriguez Espinosa, 2001) 
since it is dependent upon the orientation of the torus to the line of 
sight. This interpretation thus provides direct proof of the unification 
of the Scyfert classes. 

A test of unification schemes is offered via observations in hard 
X-rays that can penetrate an obscuring (Compton-thin) torus. Thus, 
in the simplest unification models, if X-ray emission is not beamed, 
the mid-infrared continuum arising directly from AGN heating of the 
torus or dust in the narrow line region is determined by the geo- 
metric distribution of the obscuring material and the aspect angle. 
After isolating the mid-infrared continuum due to the active nucleus 
alone, Lutz et al. (2004a) found no distinction between Sey Is and 2s 
with the X-ray/mid- infrared flux ratio, contrary to the predictions of 
unified schemes. The authors suggest the apparent uniformity is due 
to the dilution of anisotropic emission by a predominant mid-infrared 
component on extended scales that emits isotropicahy. 

Another result that is difficult to reconcile with the torus model 
arises from ISO SED analysis of a sample of hard X-ray selected sources 
comprising both type 1 and 2 Scyfert galaxies. Kuraszkiewicz et al. 
(2003) showed that while Seyfert 2s are redder in their optical colours 
than Seyfert Is, their mid- to far-infrared emission, as traced by the 
25 to 60/im colour, remains similar. More complex geometries for the 
obscuring material are required to explain both the difference in optical 
colours and the similarity in the infrared. Clumpy and extended models 
are suggested where a disk of only a few hundred parsecs can reproduce 
the IR SED without invoking a starburst component (see Elitzur et al., 
2003 and Sect. 6.6). 



5.6. Evolution of compact steep spectrum radio galaxies 



ISO data contributed towards understanding whether the compact- 
ness of steep-spectrum radio galaxies was due to age or evolution of 
environmental conditions that inhibits the growth of the radio lobes. 
It was expected that dense material causing this inhibition should be 
accompanied by far-infrared emitting dust. However, Fanti et al. (2000) 
found no far-infrared excess in compact steep spectrum radio galaxies, 
in comparison to a control sample of 16 extended radio galaxies, which 
favours a youthful scenario to explain their compactness. 
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Figure 3. Illustration of the unified scheme of quasars complemented by its 
manifestation in the infrared SEDs of quasars (Haas et al. 1998). 



5.7. Constraints from line profiles 
5.7.1. Probing the NLR 

The emission line profiles of Seyfert galaxies are known to exhibit blue 
asymmetries and blue-shifts with respect to the systemic velocity in 
their optical spectra. The profiles provide constraints on the NLR. 
Narrow lines detected in the mid-infrared have the advantage of be- 
ing much less sensitive to extinction than optical lines and therefore 
provide a relatively unobscured view to the NLR. The similarity in line 
shapes between the optical and the mid-infrared provides constraints 
on the nature of the NLR and models how such asymmetries were 
produced. Sturm et al. (1999) are able to rule out simple radial flows 
plus dust scenarios and they propose that a central geometrically thin, 
but optically thick, obscuring screen located close to the nucleus is 
present in NGC 4151 to explain the observed asymmetry. The mid- 
infrared line profiles of NGC 1068 are due to a combination of an 
highly ionised outflow and extended emission of lower excitation, and 
remaining blue shifts attributed to an intrinsic asymmetry in the NLR 
or very high column density in the line of sight (Lutz et al., 2000b). 
For a larger sample, Sturm et al. (2002) show that the differences 
in line profiles depend upon the ionisation potential. Circinus shows 
symmetric low-ionisation and asymmetric high-ionisation lines while 
the spectrum of NGC7582 shows the reverse. These differences suggest 
that the NLR properties affecting individual galaxies are particular to 
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that galaxy. Differences including extinction, asymmetric distributions 
of NLR clouds, orientation must all influence the line shapes (Sturm 
et al., 2002). 

5.7.2. Probing the BLR 

Lutz et al. (2000a) searched for broad components in the H-rccombination 
lines of NGC 1068, since the detection would imply that at least part 
of the BLR is accessible in the mid-infrared (looking through the ob- 
scuring putative torus). The lack of detection places a lower limit in 
both extinction (Ay > 50mag) and hydrogen column density (Nh > 
10^^ cm~^). Although little can be stated regarding favoured torus mod- 
els or constraining gas properties based upon this information. Bra 
was identified as being the preferred indication for future studies and 
Lutz et al. (2002) found broad line components in approximately 1/4 
of their sample of Seyfert galaxies from ground based high resolution 
spectroscopy. 



6. Ambiguous Sources: Ultra- & hyperluminous infrared 

galaxies 

As we have previously discussed, ULIG samples contain galaxies of both 
starburst and AGN type. Establishing whether these types are linked 
in evolution or nature is a key issue in the study of active galaxies 
raising several questions: what is the role of star formation and black 
holes; which is/are present; do super-massive black holes and starburst 
activity concurrently exist in all ULIGs and how much do they con- 
tribute to the bolometric luminosity? Prior to ISO, few observational 
tools were available to probe through the dust obscuration to reach 
the source of the infrared power. ISO's spectroscopic and photometric 
capabilities enabled significant advances in this field, largely elucidating 
the previously termed 'starburst- AGN controversy' which is discussed 
in this section. 

6.1. Quantitative Spectroscopy 

6.1.1. Fine Structure Line Diagnostics 

Fine structure lines directly trace the excitation states of the ISM from 
which they originate. The radiation fields generated by AGN/QSO 
greatly exceed that of starbursts and thus can excite interstellar gas 
to higher ionisation species (as described in Sect. 2.2.2). The difference 
in excitation between AGN and starbursts is the basis of fine structure 
line diagnostic diagrams. Such excitation diagrams are well established 
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Figure 4- Diagnostic diagrams based upon fine structure lines from Sturm et al. 
(2002). Fig. a: [NcVI]/[OIV] vs [NeVI]/[NeII]: Diamonds - Soy Is, Asterisks - Scy 
2s: the composite Seyfert galaxies that exhibit UIBs are encircled and are mostly 
separated from the remaining galaxies. Fig. b: [OIV]26/um/Br/3 vs [Si II]34/im/Br/3: 
Diamonds - AGN asterisks - Starbursts: the starbursts occupy a different regions of 
the diagnostic diagram than the AGN. 



in the optical (Veilleux and Osterbrock, 1987; Kewley et al., 2001), for 
which Sturm et al. (2002) constructed the first infrared analogues (Fig. 
4) that successfully distinguish starbursts from AGN. 

Diagnostic diagrams from Sturm et al. (2002) are shown in Fig. 4 
both involving the fine structure line [OIV]25.91jUm. While this line 
was expected in the high ionisation potential of AGN, it was also 
often weakly detected in star-forming galaxies and ULIGs (Lutz et al., 
1998a). In our Galaxy, the [OIV]25.91//m line has not been observed 
in HII regions surrounding young, hot, massive stars. However, its 
detection in ISO spectra of starbursts and the fact that it was spatially 
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resolved in the low excitation starburst M82 implies it is plausibly 
produced by ionising shocks or extremely hot ionising stars (Lutz et al., 
1998a; Schaerer and Stasihska, 1999). Its diagnostic capability thus 
arises from the fact that the measured strength of the line in starburst 
galaxy spectra is at least ten times fainter than those measured in spec- 
tra of AGN. This absolute difference in the strength of [OIV] relative to 
a low excitation starburst line such as [Nell] provides a straightforward 
indicator of AGN versus starburst activity (Sturm et al., 2002). 

Combining the high:high excitation fine structure line ratio [NeVI]/[OIV] 
with the highdow excitation fine structure line ratio [Ne VI] / [Nell] , 
provides a basic diagnostic of composite sources (Fig. 4a). By using 
the limiting values of the highdow excitation ratio [OIV] /[Nell] (always 
<0.01 in starbursts and 0.1<AGN<1.0) Sturm et al. (2002) were able 
to estimate the fractional contribution of an AGN to the bolonictric 
luminosity through a simple mixing model. Fig. 4b shows how the high 
and low ionisation indicators normalised to the H-recombination line 
Br/3 (an indicator of the Lyman continuum rate) provides an analogue 
to the Veilleux and Osterbrock (1987) diagnostics that are insensitive 
to obscuration and will be extremely useful for classifying ULIGs ob- 
served with future infrared spectrographs. An equivalent diagnostic 
at longer wavelengths can be constructed using [CII]158/im/[OI]63/im 
versus [OIII]88//m/[OI]63//m that can separate starburst, AGN and 
PDR contributions (Spinoglio et al., 2003). 

As the spectroscopic sensitivity of ISO permitted the detection of 
H-recombination lines or useful limits on the [OIV] emission in only 
the brightest sources, alternative tools using more commonly detected 
spectral features (e.g. UIBs and bright fine structure lines) were used 
during earlier ISO science analysis. Genzel et al. (1998) combined exci- 
tation and UIB strength to asses the nature of a large sample of AGN, 
starbursts and ULIGs. Laurent et al. (2000) constructed a diagnostic 
based upon the mid-infrared continuum shape to UIB strength (Sect. 
6.1.2). However the use of UIBs is limited and presents problems in 
interpretation (Sect. 6.1.3). 

6.1.2. UIB-based diagnostics 

The predominance of UIBs in starbursts (and their deficiency in AGN, 
or low metallicity environments) naturally lends the ratio of UIB fea- 
ture flux (normalised to the mid-infrared continuum) to be a suit- 
able tracer of star formation activity (Verstraete et al., 1996; Forster 
Schreiber et al., 2004; Peeters et al., 2004). For normal and starburst 
galaxies the UIB to continuum ratio is seen to be high in the nuclear 
star- forming regions and to decrease with radial distance (Mattila et al., 
1999; Dale et al., 2000), indicating a weakening ISRF. Conversely, a 
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strong enhancement in the continuum to hne ratio is seen in the vicinity 
of AGN, where the continuum is strong and UIBs weak as the features 
suffer dilution due to the absolute high level of the continuum and/or 
destruction of their carriers (Lutz et al., 1998b). The effectiveness of 
this indicator is clearly seen in Centaurus A, where spatially resolved 
ISOCAM-CYF observations of the AGN-dominated nucleus and the 
extended starburst region (Mirabel et al., 1999) show striking differ- 
ences in their mid-infrared emission. Together with the absence of UIB 
emission features, the spectral index of the continuum is steeper in AGN 
than in the surrounding circumnuclear star forming regions. These 
spectral differences enable, at least in principle, discrimination between 
the emission characteristics of the two distinct physical mechanisms. 

Genzel et al. (1998) demonstrated the effectiveness of the 7.7/im 
UIB-to-continuum ratio versus a tracer of hardness (such as the ratio 
of the high excitation potential fine structure line [OIV] to the low 
excitation line [Nell] (or [SIII])) to classify ambiguous sources. This 
'Genzel-diagram' (Fig. 5a) enabled the classification of ULIGs for which 
the underlying fuelling mechanisms were unclear in the pre-ISO era. 
The separation between AGN, starbursts and the intervening ULIGs is 
unequivocally depicted. The well known low-excitation starburst M82 
is far separated from the AGN NGC 4151. These extremities more or 
less define an appropriate mixing line that begins with 0% AGN in the 
lower right corner to 100% in the upper right. Thus placing a ULIG 
on this plane provides an indication of its AGN fraction. The line of 
equivalent mixing is marked on the diagram, 4 of the 13 ULIGs shown 
fall to the left of that line implying a bolometrically dominant AGN is 
needed to explain their emission. 

Similarly, using the ratios of 6.2fim UIB-to-continuum vs. warm-to- 
hot continuum, Laurent et al. (2000) constructed a diagnostic diagram 
that can quantitatively estimate the fractional contribution from AGN, 
H II regions and PDRs. Laurent et al. (2000) showed how a mid-infrared 
spectrum may be quantitatively decomposed into these three compo- 
nents through spectral template fitting (see Fig. 5b), using the pure 
HII spectrum of M17 to fit the strong VSG continuum in starbursts 
(Cesarsky et al., 1996b), an isolated PDR spectrum of the reflection 
nebula NGC 7023 (Cesarsky et al., 1996a) and the nuclear spectrum 
of Centaurus A representing the hot AGN spectrum (Laurent et al., 
2000). Such diagnostics are extremely valuable for assessing the fu- 
elling mechanisms of ambiguous sources such as ULIGs (Laurent et al., 
2000; Tran et al., 2001) and can reveal the presence of obscured AGN in 
the mid-infrared. This diagnostic tool is most successful when used with 
spectra of high signal-to-noise and with sufficient coverage of the mid- 
infrared range [i.e. to best determine the continua shortward (<5//m) 



Page 38 



Obscured "Activity": AGN, Quasars, Starbursts and ULIGs 



39 




d I I I I I I I I I z 

ft.l I 10 



relative strength of 7.7[jm PAH fca^turc 

Fig a 



1 — ' — ' — ' — I — ' — ' — ' — I — ' — ' — ' — r 




D.OOOl I I . I I I I I I I I I I I I I . I I I I I i_ 

0,0001 0,0010 0,0100 

6.2^m -PAH/FIR 

Fig c 

Figure 5. UIB diagnostics from Genzel et al. (1998) (a) and Laurent et al. (2000) 
(b). Both of these diagrams show the mixing levels of starbursts and AGN in ac- 
tive sources. Finally the MIR/FIR diagnostic of Feelers et al. 2004 is shown for 
comparison (c). 
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and longward (>12;um) of the family of strong UIB features] that en- 
ables accurate decomposition of the H II PDR and AGN templates and 
reduces problems associated with blended mid-infrared absorption and 
emission lines (Spoon et al., 2002). A follow up of this method, on a 
larger sample of galaxies was presented by Peeters et al. (2004) with 
similar results. The authors also construct a new mid- and far-infrared 
combined diagnostic that is useful for differentiating obscured sources 
(Fig. 5c). In integrated galaxy spectra, where no spatial information is 
available, e.g. for compact and/or distant sources, decomposition of the 
spectral characteristics is the only means to identify the power sources 
constituting the mid-infrared emission (e.g. Vigroux, 1997, Rigopoulou 
et al., 1999; Lutz et al., 1998b; Mouri et al., 1998; Genzel et al., 1998; 
Laurent et al., 2000). 

6.1.3. Limitations 

The use of the UIB strength as a diagnostic has several limitations. 
One originates from the fact that measuring the strength of the l.lfim 
UIB feature in the wavelength range of ISOPHOT-S (2-11/xm) was 
challenging as ULIGs also display strong silicate absorption at 9.7um 
and tracing the continuum was often extinction dependent. Another 
limitation of relying on UIB strength originates from its sensitivity to 
the radiation environment. Its absence in hard radiation fields in the 
vicinity of massive stars or AGN weakens its effectiveness as a tracer 
of star formation. Finally, as discussed by Laurent et al. (1999, 2000), 
ISO's limited spatial resolution in the mid-infrared, may introduce a 
scatter into the diagnostic diagram, especially for distant sources with 
angular size smaller than the instrument beams. This results in a dilu- 
tion of the AGN signatures as a larger fraction of the circummiclear star 
forming regions of the host galaxy disk enter into the beam. For these 
reasons, while UIB diagnostics were relatively successful in particular 
for faint sources. The diagnostics formed purely from fine structure 
lines, which directly trace differing excitation regions within galaxies, 
offer far cleaner and more effective means of source-type segregation 
and should be preferred for classification of faint and high-redshift 
galaxies that will be observed by future telescopes. 

6.1.4. Broad-band colour diagnostics 

The well known diagnostics based on the broad band infrared colours 

determined from IRAS data (e.g. the IRAS 60/im/25/im colour repre- 
sented the warm/cold emission ratio) have been extended to include 
ISO photometry (Spinoglio et al., 2002) and show the separation of 
'warm' AGN (Sect. 5.2), 'cold' starburst (Sect. 4) and 'very cold' (Sect. 
3) dust dominated systems (see their Fig. 14). While in the far-infrared. 



Page 40 



Obscured "Activity": AGN, Quasars, Starbursts and ULIGs 



41 



the SEDs of Sey Is and 2s are indistinguishable, the diagram of Spinogho 
ct al. (2002) is based on the fact that the SEDs of Seyls are flatter in 
the mid- to far-infrared than Scy 2s. Using averaged spectra from this 
sample, the authors were able to produce redshift tracks on colour- 
colour diagrams for the Spitzer MIPS and IRAC filters as a means of 
determining infrared photometric redshifts. Even though these tracks 
were calculated without evolutionary considerations, the diagrams are 
useful tools to differentiate AGN (Seyfert 1 &2s, QSOs) dominated 
sources from starbursts. 

6.2. The nature of ULIGs 

Using these techniques, the majority of ULIGs arc found to be predom- 
inantly powered by starbursts. For large samples of ULIGs Lutz et al. 
(1998b), Rigopoulou et al. (1999) and Tran et al. (2001) corroborate 
the early result of (Genzel et al., 1998) and find that ~75-85% of their 
sample do not contain energetically dominant AGN. Starburst activity 
can explain the bulk of the bolometric luminosity although the presence 
of a weak or heavily obscured AGN cannot be excluded. While the 
fraction of AGN dominated sources is 15-25%, the fraction of AGN 
hosting ULIGs is 20-25%. The AGN dominated sources have 'warm' 
IRAS colours {825/8^0 ^ 0.2) confirming previous IRAS results on 
ULIGs with optical or near-infrared spectroscopic classification. More- 
over, a comparison of the mid-infrared classification with good quality 
optical spectroscopy is remarkably consistent providing both optically 
classified H II regions and LINERs are both ascribed to the mid-infrared 
starburst class (Lutz et al., 1999). 

6.3. The nature of LINERs 

The nuclei of a large fraction (as much as one-third) of local galaxies 
display spectra with weak emission lines of lower ionisation than is 
seen in typical AGN, e.g. Seyfert nuclei. Furthermore, the low ionisation 
species present can not be explained through photoionisation by normal 
stars. The population of low ionisation nuclear emission-line regions 
- LINERs (Heckman, 1980) is well known to be mixed. Because of 
their prevalence, establishing the fraction of AGN or star-forming like 
LINERs has strong implications on our understanding of galaxies. ISO 
spectroscopic results clearly indicate that, at least in the infrared, their 
properties are more similar to starbursts (Lutz ct al., 1999). In this 
case, the low ionisation emission lines are ascribed to ionising shocks 
related to starburst driven winds (Lutz et al., 1999; Sugai and Malkan, 
2000). Interestingly though, in an infrared study of LINERs, Satyapal 
et al. (2004) found as much as two-thirds of the sample with hard 
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X-ray counterparts have compact X-ray morphologies consistent with 
those found for AGN. The remainder with scattered X-ray emission 
are brighter in the infrared than the compact ones. The authors found 
that most LINERs displaying high excitation mid-infrared lines also 
show compact X-ray morphology, even though a few systems deviate 
from this (i.e. NGC404 and NGC6240). These studies indicate that 
infrared bright LINERs are most likely members of a shock-dominated 
population whereas optical selection may trace the low luminosity AGN 
component. The fact that LINERs occupy a locus in ionisation di- 
agrams situated between classical starbursts and AGN suggests an 
evolutionary connection whereby LINERs may represent a phase when 
the AGN is deeply buried and the unrelated starburst dominates the 
detected emission (Lutz et al., 1999). 

6.4. ULIG SEDs 

An alternative classification method arises from an analysis of ULIG 
SEDs that include ISO data. This includes fainter and/or more lumi- 
nous ULIGs and HyLIGs that were inaccessible to ISO's spectrometers. 
Often in conjunction with multi-wavelength data from other telescopes, 
ISO data has enabled classification of ULIGs through decomposition of 
their SEDs. A statistical study of 20 ULIGs (Klaas et al., 2001) clearly 
showed a progression of SED shapes. These were broadly divided into 
two types. The first (type A) show a power law increase from 1 to 
50fim and are classified as Seyferts from optical spectroscopy. They 
also have warm infrared colours {hbum / feo/xm > 0.2). Type B sources 
display a relatively flat SED from 1 to 10/im, that is followed by a 
stcpp rise towards the peak of emission in the far-infrared. These SED 
characteristics are displayed by ULIGs with optical spectroscopic clas- 
sifications of Seyfert 2, LINER, Hll/starburst. These are 'cool' ULIGs 
with infrared colour f25^m / feo/tm < 0.2. The grouping of LINERs in the 
'cool'-starburst-like ULIG class is consistent with the analysis of mid- 
infrared versus optical spectroscopic classification performed by Lutz 
et al. (1999) discussed in the previous section. This, together with the 
near-infrared colours (see Figure 6 in Klaas et al. (2001)) of ULIGs, sug- 
gest that the far-infrared emission from LINER-ULIGs is of starburst 
origin. Klaas et al. (2001) interpreted the emission of 'warm'-ULIGs 
(type A) to arise from dust heated to high temperatures directly from 
the central AGN, which in turn dominates the mid-infrared emission. 
The second SED feature is ascribed due to thermal re-radiation by 
dust heated by starburst activity. The indistinguishable far-infrared 
SEDs of AGN- and starburst-dominated ULIGs suggests that the far- 
infrared emission largely comes from less active or shielded regions that 
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are generally not heated by the AGN (Klaas et al., 2001). Thus Klaas 
ct al. (2001) proposed a three-stage dust model consisting of hot AGN 
heated dust, warm starburst heated dust (50K>T>30K) and cold, pre- 
starburst or cirrus like dust (30K>T>10K). The general consensus 
from ISO spectroscopic results of warm ULIGs being AGN-dominated 
and cold ULIGs starburst-dominated is corroborated by the findings of 
Klaas et al. (2001). 

6.5. HyLIGs and the increasing role of AGN at high 

LUMINOSITY 

The increasing importance of AGN with increasing luminosity (Shier 
et al., 1996; Veilleux et al., 1997) has been confirmed by ISO obser- 
vations of ULIGs (Lutz et al., 1998b; Tran et al., 2001; Charman- 
daris et al., 2002). Tran et al. (2001) quantified that the change from 
starburst- to AGN-dominated systems occurs at logio(L/K)~12.4-12.5L0 
While ISO results show that lO^^L© can be powered by star forma- 
tion alone (Genzel et al., 1998; Lutz et al., 1998b; Rigopoulou et al., 
1999; Tran et al., 2001), beyond IO^^Lq an AGN contribution seems 
inevitable. Moreover, the production of such a luminosity would re- 
quire star formation rates ^ fewxlO'^MQyr""^. Such rates require huge 
concentrations of molecular gas to be present (most likely the result of 
mergers) and highly ineffective negative feedback (Takagi et al., 2003). 

The prevalence of AGN in HyLIG samples partially arises from a 
selection eS'ect. These rare luminous sources have been discovered with 
heterogeneous selection methods, most commonly through a correlation 
of known quasar or radio catalogues with IRAS or ISO surveys resulting 
in the known population being biased towards AGN. Constraining the 
properties of this sample and full consideration of the biases involved 
is unfortunately not possible due to the paucity of known HyLIGs. 
Nevertheless ISO observed several hyperluminous galaxies in a range 
of programs. 

Only a few HyLIGs were sufficiently bright to permit low resolution 
spectroscopic observations. Analyses of IRAS 09104-1-4109 (Tran et al., 
2001) and IRAS F15307-F3252 (Aussel et al., 1998) showed that the 
mid-infrared emission predominantly originates from very compact re- 
gions <100pc (Taniguchi et al., 1997) with lower limits on the size of 
the emission region of 6pc for F15307 (Aussel et al., 1998) and 5/-y/m 
(m=magnification) for the Cloverleaf. These quasars show hot gas at 
temperatures in excess of 400K. They exhibit low UIB-to-continuum 
ratios expected for AGN. For F15307 SED modelling showed that a 
starburst component is required to fit this source fuelling 45% of its 
luminosity, although this poses a problem for the non-detection of CO 
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(Verma et al., 2002). In a sample of four HyLIGs including HyLIGs 
without clear AGN features, a SED analysis indicated an AGN was 
required to explain the IRAS-ISO SEDs (Verma et al., 2002). One of the 
highest redshift sources observed by ISO was the powerful radio galaxy 
8C 1435+635 (z=4.25). It was shown to be hyper luminous and massive 
but the infrared-mm emission was solely ascribed to star formation and 
extremely high star formation rates were implied (Ivison et al., 1998). 

Several HyLIGs are also well known PG quasars that do not conform 
to being dust free naked quasars (Haas et al. , 2000b) . Haas et al. (2003) 
show a schematic view of the evolution of SEDs encompassing various 
types. A significant AGN contribution is seen but PG quasars lack 
the cooler extended dust component seen in lower luminosity sources. 
Patchy or clumpy disk torus models may explain all of the far-infrared 
emission, without invoking starburst contributions. 

6.6. Starburst components in AGN-hosting galaxies 

Is concurrent star formation a requisite to explain the SEDs of QSOs 
displaying a far-infrared excess? Even the cold far-infrared-sub-mm 
emission of AGN dominated sources may be fuelled by the central 
obscured active nuclei itself. While ISO data cannot spatially resolve 
such emission, invoking strong star formation is not always necessary 
to describe the SEDs of AGN and ultraluminous quasars (see below) . A 
clumpy distribution of obscuring material around the black hole allows 
high energy photons to escape and heat extended dust (Nenkova et al. 
(2002), Elitzur et al. (2003)). Alternatively warped, thick and extended 
disks (on parsec scales) have been proposed to explain the emission (e.g. 
Granato and Danese (1994)). 

There is clear evidence that some AGN do not host bolometrically 
important starbursts. For example, the archetypal Seyfert 2 galaxy 
NGC5252, shows no signs for starburst activity and the bulk of its 
far-infrared SED can be explained by the heating of the cold dust by 
the general ISRF while a small fraction is ascribed to a very large 
disk (kpc scale) surrounding the active nucleus (Prieto and Acosta- 
Pulido, 2003). The infrared emission of narrow- line Seyfert 1 galaxies 
measured by ISOPHOT can be accounted for by re-radiation of the 
soft X-ray component (Polletta and Courvoisier, 1999a; Polletta and 
Courvoisier, 1999b). Similarly, at higher AGN power, starburst activity 
is not a necessary requisite to explain the full SEDs of a large sample of 
3CR sources (Siebenmorgen et al., 2004a) or PG quasars (Haas et al., 
2000b; Haas et al., 2003) for which AGN heating plus extended cold 
dust suffices (also see Rocca-Volmerange and Remazeilles, 2005). In a 
sample of 22 quasars, starburst emission does contribute to the far- 
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infrared emission but to a level no more than 27% (Polletta et al., 
2000). 

Thus while concurrent starbursts may be present, ISO SEDs of op- 
tically classified AGN and quasars indicate that they are typically not 
bolometrically dominant. The detection of spatially resolved star burst 
tracers in the host galaxies of quasars or AGN-like ULIGs and HyLIGs 
would clarify this issue. 

6.7. AGN COMPONENTS IN STARBURST-DOMINATED GALAXIES 

Infrared diagnostics were largely successful in assessing the nature of 
most infrared sources. However, how strong a limit can one place upon 
the absence of an AGN in the mid-infrared? At short mid-infrared 
wavelengths < 10/xm the strong VSG continuum related to the star- 
formation in starburst dominated sources will overshadow the weaker 
AGN continuum, particularly in Sey 2s where the continuum due to the 
AGN will suffer high obscuration (Laurent et al., 2000; Schulz et al., 
2000). This results in an AGN being identified through the Laurent 
et al. (2000) diagnostic only if it dominates over star formation. Galax- 
ies such as IRAS 00183-7111 and Arp 220 are highly obscured even in 
the mid-infrared (Sturm et al., 1996; Tran et al., 2001; Spoon et al., 
2004) , Arp 220 may be optically thick even in the far-infrared (Fischer 
et al, 1997; Haas et al., 2001; Gonzalez- Alfonso et al., 2004)). This 
substantial obscuration by dust may cause any AGN to be completely 
hidden (Schulz et al. , 2000) , a fact which is compounded by the possibil- 
ity that AGN tracers could in some cases be completely overshadowed 
by circumnuclear activity (Laurent et al., 1999; Dennefeld et al., 2003). 
Without spatial resolution, the possibility of detecting extremely ob- 
scured or weak AGN in, at least the mid-infrared, appears to be low 
(Laurent et al., 2000; Schulz et al., 2000). 

6.8. Newly detected ULIGs 

As well as the ULIGs and HyLIGs described here which were known 
prior to the launch of ISO, several were discovered primarily in ISO 
surveys. Among the many ULIGs detected in the deep ISO surveys, 
a hyperluminous QSO (Morel et al., 2001) was discovered in ELAIS, 
and a highly extincted (Ay ~ 1000) gas and dust rich elliptical was 
detected in the /S'OPiJOT Serendipitous Survey (Krause et al., 2003). 
One might expect the population of heavily obscured extremely red 
objects (EROs) to be luminous in the infrared. Indeed, ultraluminous 
emission (and high star formation rates) was determined from the ISO 
SEDs of the ERO HRIO with a SED similar to that of Arp220 (El- 
baz et al., 2002) and a sub-mJy radio ERO from the Phoenix survey 
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(where coeval star formation and AGN are probable (Afonso et al., 
2001)) through targeted observations. Tentative infrared detections of 
a gamma ray burst were also reported to be at ULIG luminosities 
(Hanlon et al., 1999; Hanlon et al., 2000) originating from either the 
host galaxy or the GRB infrared afterglow. 

6.9. Evolution of ULIGs 

The basic similarity in luminosity and space density of ULIGs and 

quasars, and the fact that unbiased samples of ULIGs contained both 
star forming and AGN members, led to the suggestions that the two 
populations are linked through evolution (Sanders et al. (1988), also 
see Sanders and Mirabel (1996) and references therein). The example 
of the composite galaxy NGC 985 is used to discuss possible ULIG 
formation and destruction scenarios ( Appleton et al. , 2002) . Some pop- 
ular evolutionary models include the following that were recently well 
summarised by Li'pari et al. (2003): 

(1) merger - giant shocks - super-starbursts + galactic winds - elliptical 
galaxies (Joseph, 1999) 

(2) merger - H2 inflow (starbursts) - cold ULIRGs - warm ULIRGs+QSOs 
- optical QSOs (Sanders, 1999) 

(3) merger/s - extreme starburst + galactic wind (inflow + outflow) - 
Fell/BAL composite IR-QSOs - standard QSOs and ellipticals (Li'pari 
et al., 2003) 



6.9.1. Merger Stage and Luminosity 

Merging activity that initiates associated starbursts leading to the 
ultraluminous behaviour in the early stages is common to all of the sce- 
narios above. Thus the merger stage should be related to the luminosity 
evolution of ULIGs. Furthermore, if quasars are the end product of 
ULIG evolution and the merger stage indicates an evolutionary status, 
then more advanced mergers should appear more AGN-like. 

The near- to mid-infrared energy distributions of ULIGs are dom- 
inated by compact nuclear regions (Charmandaris et al., 2002). This 
confirms results from high resolution ground-based observations (Soifer 
et al., 2000) that show 30-100% of the mid- infrared flux can be ac- 
counted for by the central 100-300pc. Based on ISO imaging of a 
merging sequence of active galaxies, including starbursts, LIGs and 
ULIGs, Charmandaris et al. (2001a) suggest an evolutionary sequence 
that is traced by the 15/im/7;um colour that monotonically increases 
with star formation activity, from 1 (quiescent) to 5 (extreme) , as galax- 
ies evolve from pre-starburst to the merging starburst phase. The ratio 
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decreases again to 1 for the post starburst phase in evolved merger 
remnants [cf local less luminous post-starburst colours of NGC 7714 
(O'Halloran et al., 2000) and NGC 1741 (O'Halloran et al., 2002)] The 
fact that the extra-nuclear activity typically contributes only a minor 
fraction to the bolometric luminosity of ULIGs (in comparison to LIGs 
(Hwang ct al., 1999)) implies that simple scaling in mass and luminosity 
does not link the two populations, but evolutionary differences (such 
as merger stage, or presence of an AGN) could play a role. It is unclear 
whether this hardening of the mid-infrared continuum is due to an 
AGN contribution. The amount of molecular gas may also affect the 
spectral feature differences seen between the populations, such as the 
presence of ice features (Spoon et al., 2002). Increasing compactness 
of sources with increasing merger stage shown by Hwang et al. (1999) 
implies that compact ULIGs are in a more advanced merger state than 
LIGs. 

However, based upon a comparison of the NIR morphologies of merg- 
ing ULIGs and infrared spectroscopic tracers, Lutz ct al. (1998b) and 
Rigopoulou et al. (1999) do not find any relation of increasing AGN- 
like activity with merger stage. In addition ULIGs appear to not show 
a trend of increasing luminosity with merger stage again contrary to 
expectations of these evolutionary scenarios (Rigopoulou et al., 1999). 
Indeed Rigopoulou et al. (1999) suggest that early merger stage ULIGs 
may be exhibiting maximal luminosities (also see Murphy et al., 2001). 
Moreover, the mass of the warm, cold and very cold dust components 
found in ULIG SEDs show no relation to merger stage (Klaas et al., 
2001). The evolution or merger-stage may not be traced well simply 
by nuclear separation alone as a merger stage indicator. More complex 
differences in the mass distribution and dynamics of ULIGs more likely 
influence the luminosity and appearance of AGN. 

6.9.2. Dust evolution in quasars 

If ULIGs are related to quasars as the dust obscured members of that 
class, understanding the dust properties of quasars also serve to test the 
evolutionary connections between quasars and ULIGs. In the Sanders 
et al. (1988) scenario a quasar is preceded by a dusty ULIG-phase 
and as a result it is unlikely that the obscuring dust will suddenly 
disappear, some traces in the infrared should remain. Haas et al. (2003) 
investigated a sample of 64 PG quasars and identified a range of SED 
shapes that could mark a transition between states. The SED analysis 
showed that the central regions of PG quasars suffer low extinction 
(Ay < 0.3) with an optical slope aopt that is independent of infrared 
properties like the near- to mid-infrared slope ajR. Therefore, with 
regard to the unification schemes a nearly face-on view onto the PG 
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quasars was assumed. The observed variety of SEDs were grouped into 
physically meaningful classes, that reflect the amount and distribution 
of the reprocessing dust around the AGN according to evolution of the 
quasar (Haas et al., 2003). Their new evolutionary scenario begins with 
dissipative cloud-cloud collisions and angular momentum constraints 
that create the organisation of dust clouds into a torus/disk like con- 
figuration. Starburst activity is initiated giving rise to a far-infrared 
bump but it becomes increasingly overpowered by the central AGN 
which results in increasing mid-infrared emission and steepening of the 
mid- to far-infrared spectral slope. As the AGN continues to grow, the 
mid-infrared remains high and the mid- to far-infrared slope decreases 
until the black hole begins to starve and a a decline in the mid-infrared 
and far-infrared emission is observed. Figure 12 in Haas et al. (2003) 
depicts this scheme, the SEDs are arranged along the expectations for 
such an evolutionary scheme. 

6.9.3. Relation of broad absorption line quasars to ULIGs 
Broad absorption line quasars are a rare class of the optically selected 
quasar population (~12%) but have a higher frequency in far-infrared 
luminous, strong FeH and weak [OHI] emission line selected samples. 
Both orientation and age arguments are consistent with the origin of 
the broad absorption features with no clear resolution which is the 
most likely cause. Elvis (2000) suggest the broad absorption arises from 
viewing a single funnel-shaped outflow of absorbing material directly 
down the flow. It has also been proposed that these sources are quasars 
viewed at an orientation directly along the radial surface of the torus, 
where winds along the surface give rise to the observed absorption. 
Alternatively they may be a young phase in quasar evolution. The 
latter together with the high occurrence of BALQSO in IRAS selected 
samples suggests that they may be a transition stage between ULIGs 
and QSOs (Voit et al., 1993; Lipari, 1994; Egami et al., 1996; Canalizo 
and Stockton, 2001; Lfpari et al., 2003). In this scenario very young 
QSOs are ejecting their gaseous envelopes at very high velocity close 
to the initiation of the active phase of the black hole (Hazard et al., 
1984). 

Two BALQSOs were serendipitously discovered in ISO surveys: the 
FeLoBAL^ ISO J005645.1-273816 was discovered in an ISOCAM dis- 
tant cluster survey (Due et al., 2002). 00 37 14.3 -42 34 55.6 was 
discovered in ELAIS (Alexander et al., 2001). For the former, high 
absorption and hot dust emission are implied by the high mid-infrared 
to UV ratio which may indicate the mid-infrared is an efficient means 

^ FeLoBALs show absorption lines from many excited states of Fe II and having 
low ionisation lines and are heavily absorbed in X-rays 



Page 48 



Obscured "Activity": AGN, Quasars, Starbursts and ULIGs 49 

to detect LoBALQSOs (Due et al., 2002). The sources 00 37 14.3 -42 34 
55.6 was the only BALQSO then known to have both infrared and hard 
X-ray detections. While it has a similar hard X-ray /MIR flux ratio to 
high redshift quasars, it differs from low-redshift BALQSOs. Alexander 
et al. attribute this to absorption. As a consequence they suggest that 
absorbed BALQSOs at high redshift should be easily detected in the 
hard X-ray bands. In addition the only radio-loud BALQSO known, 
1556-1-3517, was observed to investigate the relation between radio-loud 
quasars and the BALQSO phenomenon. Mid-infrared detections im- 
plied 4M0 of dust lie along the line of sight which most likely arises from 
the molecular torus, the NLR, or a recent starburst rather than from 
the BALQSO wind itself (Clavel, 1998). ISO data confirms that even in 
this diverse, small sample, the BALQSO phenomenon and obscuration 
(i.e. infrared emission) are connected. Lfpari et al. (2003) speculate on 
an evolutionary connection between BALQSOs and ULIGs. 



6.10. Summary 

Evaluating which of the three scenarios mentioned earlier in this Sec- 
tion is most realistic will continue to be addressed by future missions. 
However, the ISO results clearly suggest that a simple Sanders et al. 
(1988) evolutionary scenario of a merger triggered starburst event with 
an evolving obscured- naked quasar scenario to be unlikely. Rather, 
it seems that the nature of the diverse infrared luminous ULIGs and 
HyLIGs has a far more complex history. Whether a starburst or an 
AGN dominates the bolometric luminosity of a ULIG is most probably 
determined by local conditions, evolution and obscuration by dust. 
Nearby ULIGs are mostly starburst dominated systems with 'cold' 
colours. This cold fraction decreases at higher luminosities where AGN 
contributions become prevalent and galaxy SEDs are 'warmer'. ULIGs 
are often linked to high-redshift (z~2) (sub-) mm emitting galaxies 
which comprise a significant fraction of the infrared background and 
contribute significantly to the star formation density at that epoch. 
If ULIGs and sub-mm galaxies are members of the same population 
then strong evolution is implied since the locally unimportant ULIGs, 
that constitue less than 1% of the local star-formation density, must 
dominate the infrared number counts at higher redshifts (Elbaz and 
Cesarsky, 2003). The strongly evolving z^l.5 LIGs discovered in ISO- 
CAM surveys (see Elbaz et al. and Oliver et al. in this volume) that 
resolve 10-60% of the far-IR-sub-mm background may provide a step 
between. 
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7. Active Galaxies as seen by ISO: a summary 



The Infrared Space Observatory has produced a wealth of information 
regarding the infrared emission properties of the fuU range of active 
galaxies. Greatly enhancing the IRAS legacy, it has provided a phys- 
ical basis upon which the next generation of infrared telescopes such 
as Spitzcr, SOFIA, ASTRO-F and Herschel will build. Most of ISO's 
power was due to the improved sensitivity and extended wavelength 
coverage. The addition of high resolution spectroscopy and develop- 
ment of new diagnostic diagrams have enabled us to probe the ongoing 
physical conditions of the ISM in active galaxies and to quantitatively 
assess the presence and bolometric contribution of star formation or 
an AGN in infrared luminous sources. Analysis of fine structure lines 
provided an excellent set of diagnostic tools that will be better pop- 
ulated and applied to a wider variety of sources in future missions 
with even higher sensitivity. The weakness of the sulphur fine structure 
lines in starburst galaxies and of the cooling line [CII]158/um in the 
most luminous active galaxies suggests that other strong lines (e.g. 
[NeII]12.8;um) should be the are preferred indicators of star forming 
objects. The finding of IRAS that galaxies with "cold" infrared colours 
are starburst dominated and "warm" are AGN dominated was broadly 
confirmed with the new ISO data. Although, it was shown that spatially 
resolved spectroscopic information is necessary to address the issue of 
the origin of the far-infrared emission in AGN and the role of concur- 
rent star formation. ISO results on AGN in general support unification 
scenarios, but suggest that a clumpy distribution of obscuring material 
around the active nucleus may be a better representation than the 
commonly assumed torus. While the majority of ULIGs are starburst 
dominated, it now appears that an increasing fraction of AGN at higher 
infrared luminosities is present. Mergers are indeed the trigger of ULIG 
activity, but the simple generic evolution of a luminous system from a 
ULIG to an optical-QSO phase appears unlikely with more local factors 
influencing the presence of an AGN. With increases in sensitivity and 
spatial resolution, it is certain that the suite of current and forthcoming 
infrared telescopes will enhance our knowledge of the obscured active 
universe by probing fainter and to higher redshifts as well as in more 
spatial detail. 
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